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ABSTRACT 

‘  t  t  I  •'  r  ' 

In  this  paper  we  present-,*  singular  perturbation  analysis  of  the 
fundamental  semiconductor  device  equations  which  form  a  system  of  three  second 
order  elliptic  differential  equations  subject  to  mixed  Neumann-Dirichlet 
boundary  conditions.  The  system  consists  of  Poisson's  equation  and  the 
continuity  equations  and  describes  potential  and  carrier  distributions  in  an 
arbitrary  semiconductor  device. 

The  singular  perturbation  parameter  is  the  minimal  Debye-length  of  the 
device  under  consideration. 

Using  matched  asymptotic  expansions  we  demonstrate  the  occurrence  of 
internal  layers  at  surfaces  across  which  the  impurity  distribution  which 
appears  as  an  inhomogeneity  of  Poisson's  equation  has  a  jump  discontinuity 
(these  surfaces  are  called  'junctions')  and  the  occurrence  of  boundary  layers 
at  semiconductor-oxide  interfaces.  '  Wer  derive  fthe  layer-equations  and  the 
reduced  problem  (charge-neutral-approximation)  and  give  existence  proofs  for 
these  problems.  They  layer  solutions  which  characterize  the  solution  of  the 
singularly  perturbed  problem  close  to  junctions  and  interfaces  resp.  are  shown 
to  decay  exponentially  Sway  from  the  junctions  and  interfaces  resp. 

We  shoirthat,  if  the  device  is  in  thermal  equilibrium,  then  the  solution 
of  the  semiconductor  problem  is  close  to  the  sum  of  the  reduced  solution  and 
the  layer  solution  assuming  that  the  singular  perturbation  parameter  is  small. 
Numerical  results  for  a  two-dimensional  diode  are  presented. 
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SIGNIFICANCE  AND  EXPLANATION 


In  this  paper  qualitative  properties  of  the  solutions  of  the  system  of  partial 
differential  equations  which  describes  potential  and  carrier  distributions  in  an 
arbitrary  semiconductor  device  are  discussed.  The  system  consists  of  Poisson's  equa¬ 
tions  for  the  potential  and  of  the  continuity  equations  for  the  carrier  densities. 

The  problem  can  be  classified  as  elliptic  (in  the  static  case).  After  appropriate 
scaling  a  small  parameter  X2,  which  is  physically  identified  as  the  square  of  the 
normed  minimal  Debyelength  of  the  device  under  consideration,  appears  as  multiplier 
of  the  Laplace-operator  in  Poisson's  equation.  Therefore,  the  system  is  singularly 
perturbed . 

We  investigate  the  asymptotic  behaviour  of  solutions  as  the  singular  perturbation 
2 

parameter  X  converges  to  zero  and  demonstrate  the  occurrence  of  internal  layers 
(thin  regions  of  fast  variations)  in  the  potential  and  the  carrier  densities.  These 
layers  occur  at  curves  in  the  device  where  n-regions  (i.e.  regions  in  which  the  doping 
profile  is  positive)  and  p-regions  (i.e.  regions  in  which  the  doping  profile  is 
negative)  meet.  In  the  analysis  we  also  prove  existence  theorems  for  the  reduced 
problem  (X2  set  to  zero)  and  for  the  layer-problems. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive  summary 
lies  with  MRC,  and  not  with  the  author  of  this  report. 


A  SINGULAR  PERTURBATION  ANALYSIS  OF  THE 
FUNDAMENTAL  SEMICONDUCTOR  DEVICE  EQUATIONS 

* 

Peter  A.  Markowich 

I.  INTRODUCTION 

An  analysis  of  the  two -d imens i ona 1  ,  static  equations  which  describe  potential  distri¬ 
bution  and  current  flow  in  a  semiconductor  device  is  presented. 

The  basis  semiconductor  device  equations  are  (see  Van  Roosbroeek  (1950)): 


( 1 . 1 )  ( a )  d  i  v  (  e  ,  7  j, )  »  q(n  - p  -  C(x,y)) 


Poisson's  equation 


(1.  l)(b)  div(Dn?n  -  unn7Ki)  *  R 


(1.  1 )  ( c  )  div(Op?P  + upp?*)  -  R 
where  the  dependent  variables  are 


(x,y)g  ::  electron  continuity  equation 
.  hole  continuity  equation 


n  :  electrostatic  potential 
n  :  electron  density 
p  :  hole  dens i ty 


is  a  bounded  domain  in  31  representing  the  device  geometry;  4  is  the  semiconductor 


permittivity  iwhich  will  be  assumed  to  be  constant  in  the  sequel);  ,  jp  are  the 
electron  and  hole  mobilities  resp 
cients  resp.  and  q  is  the  elementary  charge. 


Q„,3  are  the  electron  and  hole  diffusion  coeffi- 
n  p 


C(x,y)  is  the  doping  profile. 

that  means  C(x,y)  is  the  difference  of  the  electrically  active  concen t ra t i on  of  do¬ 
nors  and  the  electrically  active  concentration  of  acceptors.  R  is  the  recombination 
generation  rate.  In  the  sequel  we  will  neglect  recomb i na ti on -genera ti on  effects,  that 
means  we  set  R»0.  Of  course  only  solutions  with  n>0,  p>0  are  admitted 


We  assume  the  validity  of  Einstein's  relation 

0  0„ 

(1.2)  —  =  -2.  =  u T(  »  const) 

n  p 

where  U-r  is  the  thermal  voltage. 

The  electron  and  hole  current  densities  Jn  and  Jp  3re  given  by 

I  1  •  3)  (a  )  Jn  *  q( Dn7n  -  j^npi) 

( 1  •  3  )  ;  0)  Jp  *  -q( 0p7P  +  Jpp7i) . 

Numerical  values  for  the  parameters  (for  Silicon  at  roomtempera ture) are  given  in  Table  1  . 

The  elliptc  system  of  differential  equat ions ( l .  1 )  has  to  be  supplemented  by  appropriate 
bou nda ry  c ond 1 1 i ons  f o r  „,n,o  determined  by  tne  device  under  consideration.  We  assume 
that  splits  up  into  three  dispoint  parts,  namely  s.:r,vt.  and  V.'  and  >•!  are 

f*  v  <  S  0  S  1  S  OS 

open,  is  connected  and  •  :.r-  l)^,r  >0  where  the  C,  are  closed  and  connected  arcs 

1  s  “  k  =0  * 

wtii  positive  ( one-d  imen  s  i  ona  1  )  lebesguemeasure. 
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Oirichlet  boundary  conditions  for  i>, n,p  are  given  on  s £Iq  (the  are  Ohmic  contacts)  and 
zero  Neumann  boundary  condi  tions  are  prescribed  on  3S1^S  (insulating  segments ).  3SJ0J  repre¬ 
sents  a  semiconductor-oxide  interface  occuring  in  MOS-technol ogy  (see  Sze  (1981) 
for  a  survey  on  MOS-devices) .  The  oxide  is  located  in  a  bounded  domain  >  which  is 
such  that  »nn  *  ?r?os  and 

(1.3)  a*  «  0,  (x  ,y)  £  ♦ 

holds.  The  carrier  densities  n,p  only  exist  in  5.  Usua 1 1 y  3*  spl i ts  into  three  parts, 
namely  3*c  (oxide-contact)  where  a  Oirichlet  condition  for  *  is  prescribed,  3e<s 
(insulating  segment )where a  zero  Neumann  condition  for  *  holds, and  the  interface 
3 *  has  to  be  continous  accross  3si„«.  and 

OS  OS 


(1.4) 


t  e  r  ]  .  0, 

3nos 


«( * .y) 


(x,y)€« 

(x,y)6n 


€  t 

holds  where  c  is  the  oxide  permittivity  (  —  *i)  (iflr  denotes  the  jump  of  the 

0  £  «  J  - 

function  f  accross  the  curve  r).  n  is  the  exterior  unit  normal  vector  of  3fl.  The 
condition  (1.4)  represents  the  continuity  of  the  electrical  displacement  accross 
the  semiconductor  oxide  interface. 


The  electron  and  hole  current  density  components  Jn'  >t  and  Jp'*1  (perpendicular  to 
3.1  )  vanish  on  3S10$.  This  gives  boundary  conditions  for  n  and  p  (by  using  (1.3))  at 
the  interface. 

The  Oirichlet  boundary  conditions  for  n  and  p  at  the  Ohmic  contacts  are  given  by  the 
vanishing-space-charge  condition 

(1.5)  (n  - p  -  C(x.y))|  »  0 

and  the  thermal  equilibrium  condition 

(l.e)  npi3(Jc  .  n2 

where  n^  is  the  intrinsic  number  of  the  semiconductor. 

For  the  following  we  assume  that  fl  splits  up  into  N+l  connected  subdomains  0., 

N  _  1 

('  •  'J  .’•.),  such  that  C  does  not  change  sign  in  each  of  the  0.  and  C  has  jump  - 
i  =o  __  1 

discontinuities  accross  the  curves  r.  =  (abrupt  doping),  f.nfj  *  {}  holds 

for  i  *j.  0^  is  called  a  n  region  if  C I  g  >0  and  it  is  called  a  p  region  if  C !  ^  <0. 
r  1  is  a  pn junction  if  it  is  the  joint  boundary  of  a  p  and  an n  region  and  it  is  in 
r.n  (pp)  junction  if  it  is  the  joint  boundary  of  two  n  (p)  regions. 


We  also  assume  that  the  Ohmic  contactsCj,  have  positive  distance  from  the  junctions 
ri  • 

The  performance  by  the  device  under  consideration  is  mainly  determined  by  the 
location  of  the  subdomains,  of  the  oxide  (for  MOS-devices)  and  by  the  location  of  the 


Ohmic  contacts.  The  boundary  conditions  for  th«  potential  (at  the  Ohmic  contacts) 
are 


(1.7)  (a)  *  Ufl n  ^ilck+Uk*  1f  Cke51  40,1  C'.^>0 


In  n-regions  and 


(1.7)  (b)  *ic^  *  UTm  ~lc*uk  if  Ck«a1  and  C|^<0, 

In  p-regions  where  fik  represents  the  potential  applied  to  the  Ohmic  contact  C^.  He 
remark  that  there  are  devices  which  are  such  that  not  every  n  or  p  region  "has"  an 
Ohmic  contact  (for  example  thyr 1 stors .see  Sze  (1981)). 


An  externally  appl led  potential  Uj,  Is  given  at  the  oxide  (gate)  contact 


(1.71(0  »:„c  •  u5  -  uF 

where(the  flat  band  voltage)  (J-  Is  a  constant  which  depends  on  the  semiconductor, 
on  the  oxide  contact  and  on  the  doping.  The  applied  potentials  Uk,  Uj.  are  constants, 
too . 


As  Illustration  for  the  notation  of  the  device  geometry  we  show  a  typical  MOS-transI stor 
In  figure  l. 


Figure  1.  Mos-Translstor 

There  are  two  n-regions  and  one  p-region  ( n-channe Vtransl stor) ,  three  Ohmic  contacts 
(C^ :  source-contact ,  Cg'-drain  contact,  C0  bulk  contact)  and  one  oxide  contact 
(a*j-:  gate  contact) . 

The  vertical  boundaries  of  ft  and  a  are  insulating  (sii<s  and  resp.). 

We  remark  that  the  rectangular  shape  of  a  and  t  as  shown  in  figure  1  is  a 
simplification  commonly  used  for  numerical  simulation  (see  Sel berherr{ 1980) ) .  The 
following  theory  however  is  not  restricted  to  particular  shapes  of  domains. 
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The  problem  (1.1)  can  be  put  into  a  simpler  form  by  the  transformation  (called 
Boltzmann  Statistics) 


(1.8) 


* 

Ur 


n^  e 


$  $ 

Here  u  ■  exp  (-  .  v  •  exp  (^  )  where  »n,  ♦ 

mllevels  resp.  (uT>0,  v>0  has  lo  hold).  Then 
(1.2)  ,  assuming  R»0  and  *,  to  be  constant) 

St  ^t 

( 1 . 9)  (a)  c  SA*  »  q^e  ii  -  n^e  'v-C(x.y)) 

ST 

(1.9) (b)  div  (une  7u)  •  0 

■K 

( 1 . 9)  (c)  div  (Upe  7v )  «  0. 


p  are  the  electron  and  hole  quasifer 
(1.1)  takes  the  form  (by  using 


1 

•  (x.y)  e  a. 


The  continuity  equations  (1.9)(b),(c)  are  in  self-adjoint  form. 


There  have  been  many  analytical  and  numerical  investigations  of  (1.1)  ((1.9)). 

Mock  (1972)  showed  the  existenceof  a  solution  of  (1.9)  subject  to  the  mixed  set 
of  boundary  conditions  and  he  proved  that  this  solution  is  unique  if  the  applied 
potentials  are  sufficiently  small.  He  only  assumed  C€l**(3).  Continuous  depen¬ 
dence  of  the  solutions  on  the  boundary  data  was  also  shown  in  this  paper. 

A  very  similar  existence  proof  was  given  recently  by  Sank,  Jerome  and  Rose  (1982). 

The  parabolic  semiconductor  problem  (with  homogenous  Neumann  boundary  conditions 
on  an)  was  investigated  by  Mock  (1974).  Finite  difference  methods  are  discussed  in 
Mock  (  1973),  (  1981). 


In  this  paper  we  scale  the  problem  (1.9)  appropriately  and  obtain  i  singular 
perturbation  problem.  The  singular  perturbation  parameter  \  is  the  minimal  normed 
Oebye  length  of  the  device  under  consideration. 

Using  matched  asymptotic  expansions  (as  v-0  +  )  we  demonstrate  the  occurance  of  a 
boundary  layer  in  4)  at  oxide-semiconductor  interfaces  and  the  occurance  of  internal 
layers  (in  *)  at  pn,  nn  and  pp  junctions.  u  and  v  are  the  slow  variables,  that 
means  they  do  not  exhibit  zero-order  layers. 

We  derive  the  reduced  problem  (vanishing  space  charge  approximation)  which  is  ob¬ 
tained  by  setting  the  singular  perturbation  parameter  to  zero,  and  the  (boundary 
and  internal)  layer  equations  and  give  existence  proofs  for  these  problems.  We  diss 
cuss  the asymptot i c  behavior  of  the  current  densities  Jn  ,Jp  (as  \ -0  +  )  and  show  the 
validity  of  the  asymptotic  expansions  for  the  equilibrium  problem  (zero  external 
potential  applied  at  the  contacts). 


The  singular  perturbation  approach  was  applied  to  the  one  dimensional  semiconducto 


problem  by  Vasileva  and  Stelmakh  (1977),  Vasileva  and  Butuzow  (1978); 

Markowich,  SI nghof er  ,Sel berherr  and  langer  (  1982  a,b). 

The  main  advantage  of  the  singular  perturbation  approach  is  that  it  gives 
qualitative  information  on  the  behaviour  of  the  solutions. This  a-priori  Information 
can  be  used  to  construct  appropriate  discretisation  methods  for  the  numerical  solu¬ 
tion  of  the  semiconductor  device  equations.  In  particular  efficient  mesh-strategies 
employing  only  a  reasonable  number  of  grid  points  (but  still  giving  accurate 
numerical  approximations  even  in  layer  regions)  can  be  obtained  (see  Markowich, 
Ringhofer  and  Selberherr  (1982)  and  Section  6). 

All  results  in  this  paper  also  hold  for  the  one  and  three  dimensional  static 
semiconductor  probl ems after  obvious  modification  of  assumptions  (e.g.  junctions  are 
then  represented  by  points  and  surfaces  resp.).  We  chose  the  two-dimensional 
semiconductor  problem  for  the  presentation  since  it  is  most  often  used  for  numerical 
simulation. 

The  paper  is  organized  as  follows.  In  Section  2  we  perform  the  scaling  which  leads 
to  the  singular  perturbation  problem  ,  prove  an  existence  and  a  regularity  theorem 
and  derive  a-priori  estimates  of  the  solutions.  In  Section  3  we  derive  (the  zeroth 
order  terms  of)  the  asymptotic  expansions,  in  Section  5  we  give  the  existence  proofs 
for  the  reduced  problem  and  for  the  layer  problems  and  the  equilibrium  problem  is 
discussed  in  Section  5.  Section  6  is  concerned  with  numerical  examples  and  with 
possible  extensions  of  the  theory. 

Table  1.  Numerical  Values  of  the  Parameters  for  Silicon  and  Silicon  oxide  at 
roomtemperature  T  <w  300  k. 


Parameter 

Phys leal  Meani ng 

Numerical  Value 

q 

elementary  charge 

10"l9As 

e. 

semi  conductor 

10'lZAs/Vcm 

5 

permittivity  constant 

Eo 

oxide 

permittivity  constant 

^  •  10'12As/Vcm 

el ectron  mobi 1 i ty 

103cm2/Vs 

'p 

hole  mobility 

103crn2/Vs 

electron  diffusion 

constant 

25  cfli  / s 

0. 

hole  diffusion  constant 

25  cm2/s 

p 

ni 

int  rinsic  number 

1010/cm3 

UT 

thermal  vol tage 

0.025  V 

The  numerical  values  given  for  'Jn>“pi0rl>Dp  have  to  be  understood  as  averages,  since 
these  quantities  are  generally  modelled  by  functions  of  x  and  y. 
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2.  THE  SINGULARLY  PERTURBED  PROBLEM 
We  assume  that  C  is  bounded  in  n  and  set 

(2.1)  C  »  sup|C(x.y)  |  ,  0  ■£ 

a  C 

and 

(2.2)  1  *  diam  (a). 

The  dependent  variables  are  scaled  as  follows 


and  the  independent  variables 

(2.3)  S  *  f  •  ys  *f.  (xs-ys)  eSs  u*s 

Then  (1.9)  transforms  to  (after  dropping  the  subscript  s): 

(2.4) (a)  x2d*  »  «2e*u  -  62e'*v  -  D 

{ 2 . 4 )  {  b )  div(el7u)  *  0  v  (x,y)60 

.*  i 

(  2 . 4  )  ( c  )  div(e  Vv  )  *  0 
and  (1.3)  remains  unchanged: 

(2.5)  a*  *  0,  (x,y)  €  » 

where 

(2.6)  (a)  »  (^  )  2  -  .  (b)  52  * 

holds.  xQ  is  the  minimal  Oebye  length  of  the  device. 

For  (2.4)  we  assumed  that  the  mobilities  un.up  are  constant  throughout  the  device 
(for  numerical  values  see  Table  1  ). 

The  following  theory  however  carries  over  to  the  case  that  un>up  are  smooth  and  posi¬ 
tive  functions  of  x  and  y. 

The  (scaled)  boundary  conditions  are 

( 2  .  7 )  ( a  )  ;»•<,,  =  “u-lii  *  Tv-ltj  =  t) 

3,5  is  ’'-is  iwis 

(the  (unit)  vector  J  is  perpendicular  to  tr.and  is  assumed  to  exist  almost  everywhere)  and 
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I 


(2. 7)(b)  u|c  •  e 


(2- 7)  (c)  *lc  •  In 


'*7 


v  I  (.  -  e 
**k 


D  +  /o2+4{4 


2t‘ 


«.  Uk 
; ^ 


for  k  «  0, . . ,r . 

(2.7)  Is  derived  from  (1.5),  (1.6),  (1.7)  (a.b)  by  using  (1.8). 


Boundary  conditions  on  a*  are 


(2.8)  (a)  7*.-c|  .  0 

}#1s 

(";  denotes  the  exterior  unit  normal  vector  of  ae  ) 


( 2 . 3) ( b) 

,!>.c  ’ 

UG 

UF 

07  *0 

and 

(2.3) (c) 

■  0 

(2 .8)  (d) 

[  £»7ji  •  n ) 

3nos 

*  0,  c*  *  \ 

r  11 
i  i° 

es 

(x.y)  €  0 

,(x,y)  e  ♦ 

(2.8) (e) 

7u-n|  . 

w  “  A  e 

7vni  .  0 

u  and  v  are  only  defined  In  a. 


For  modern  devices  C  »  I0*7cm*2.  With  the  realistic  value  1  »  5  x  10'^cm  and  the  numerical 
values  for  q,  es,UT  given  In  Table  1  we  get  k2  < 10'2  <<  1.  Therefore  tne  problem  (2.4), 
( 2 . 5)  ,  ( 2 . 7 )  ,  ( 2 . 3)  constitutes  a  singularly  perturbed  quasilinear  elliptic 
system  of  differential  equations  (subject  to  mixed  Neumann-01 rl chi et  boundary  and 
Interface  conditions). 

2  2-7 

The  parameter  j  <<  1,  too  (normally  4  <10  holds).  This  however  gets  compensated 

by  the  Olrlchlet  boundary  conditions  (2.7)(d),(e)  which  imply  that  S^fe^u  -e"*v)  «0(1) 

? 

at  Ohmic  contacts  as  J  -.0.  Note  that  the  potential  difference  between  an  Ohmic  con¬ 
tact  in  a  n-region  and  an  Ohmic  contact  in  a  p-region  behaves  asymptocally  (as 
iZ  -0)  like  In  -4  . 

o 

2 

We  regard  4  as  (fixed)  parameter  and  investigate  the  asymptotic  properties  of  the 
solutions  as  \  -0+.  We  will  then  show  that  the  asymptotics  are  uniform  in  5  as  long 
as  j  does  not  converge  to  zero  too  fast  (compared  to  k  ) . 
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The  scaling  factors  for  the  current  densities  Jp  are  qunE'J.j.  and  qupCu.j.  resp. 

Then  the  scaled  current  densities  are  given  by 

( 2 . 9)  (a)  Jn  -  sV?u.  (b)  Jp  •  -a2e‘*vv. 

The  scaled  carrier  densities  follow  from  (2.2) 

(2.3) (c)  n  «  S2e*u,  (d)  p  *  52e"*v. 

We  now  give  definitions  which  will  be  iieeded  in  the  sequel. 

We  denote  by  L ^ ( n )  the  space  of  q-integrable  real  valued  functions  defined  on  n  with  the 
norm 

1 

A 

"f"q,n  (/|f(*.y)|qoxdy)q 

and  by  L*,(3)  the  space  of  bounded  functions  on  a  and 

Ilf  II  .  *  sup  I  f  (x  ,y )  I  . 

*  n 

Cm ( fJ )  for  m  €  is  the  space  of  all  functions  defined  on  !i  which  together  with  their 
partial  derivations  of  order  up  to  m  are  continuous  in  n  . 

Cm(S)  is  the  space  of  all  functions  which  are  in  Cm(ii)  and  which  together  with  their 
partial  derivatives  of  order  up  to  m  are  bounded  and  uniformly  continous  in  3.  A  norm 
on  Cm(q)  is  given  by 

_  :  ,  HI 

llf|m"‘  «  max  sup  — 2 - — -  f(x.y) 

n  3X  l3y  2 

where  a*(a^!a2)  and  lal  ■a^+a2 

Spaces  of  Holder  continuous  function  Cm,s(£l),  Cm  ’ S  ( a )  for  0*3*1  and  the>r  norms  are 
as  In  Adams  (1975) . 


For  rcan  we  denote  by  H^(pur)  the  completion  of  {u  e  Cm( n) :  II u  ||  ^  p  <  »  and  u  vanishes 
in  a  neighbourhood  of  n  with  respect  to  II  .  II  m  .,  fl. 
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For  the  following  we  set  a  ■  .]  j  39  us  and  state  the  existence  theorem  which  Is 

os 

basic  for  the  theory  of  the  semi  conductor  equations. 

Theorem  2 . 1  Assume  that  0  is  defined  in  S,  0€L*(o)  and  that  0  is 

Lipschltzcontlnuous  in  a  neighbourhood  of  Also  assume  that  3A 

and  an  are  Lipschltzcontlnuous  and  piecewise  C*  and  that  the  (one* 
dimensional)  measure  of  3ftg  is  positive. 

Then  the  problem  (2.4) ,{2.5)  .(2. 7)  .(2.8)  has  a  weak  solution  (»,u,v) 
for  which  €  H*(a)  0  L°*(A)  ,  u,v  €  H^n)  fl  l“(0)  holds. 


Every  solution  ( *  ,u ,  v )  €  H*  ( A )  x  ( H*  (a)  2 ) ,  €  L°"(  A )  fulfills  the  a  priori 

estimates 

U  k  .!ii 


(  2 . 10) (a )  u  :  min 
k 


-u* 

e  <  u(x. 


y)  i  max  e 
k 

U, 


(2.l0)(b)  v  ;  min  e  *  v(x,y)  i  max  e 
k  k 

for  ( x  ,y )  6  3  and 


(2. 10)(c)  •!>_:  min  (7_ ,  inftji)  <4>(x,y)  <  max  {7*,  sup 


for  ( x  ,y )  €  a  where 


0  +  V0‘  ♦  4« \  v  ' 

(2.111(a)  In  - « -  ,  0¥  *  sup  D(x,y) 


(2. ll)(b)  *_:•  In 


0  +rOt+4S  u  v 


,0  *  i nf  D( x ,y ) 

a 


Proofs  of  slighty  weaker  existence  theorems  were  given  by  Mock  (1974)  and  Bank, 
Jerome  and  Rose  (1982).  We  therefore  only  sketch  the 

Proof :  The  a  priori  estimates  (2.10)(a)(b)  follow  immediately  by  on  application  of 
the  maximum  principle  for  H ^-sol u ti on s( see  Gilbarg  and  Trudinger  (1977), 
Chapter  8)  to  ( 2 . 4) ( b) , (c) ,  (2.7)(a),(b)  and  (2.8). 


We  set 


Q(*)  *  A2div(c*Vj,)  -  (<52e*u *  -  a2e'%*-D*) 


r  0  in  a 


m 


The  weak  formulation  of  Q(n)  *  0  subject  to  ( 2 . 7 ) ( a ) , ( 2 . 3 ) ( a ) , ( 2 . 8 ) (c ) , ( 2 . 8 ) ( c ) . 
(d)  is 


0(*,io)  •  ;(x2e*v*-?u>  ♦  («2e We*V-D*)<i>)<lxdy 
A  t 

■  /( x27*-7(p+(  42e*'u-«2e*'tlv-O)i0dxdy  +  x2  —  /7*'7ipdxdy 

Q  S  * 

■  0 

for  all  <o  €  H*  (  a  u  3!^  u  3*  ^  )  since  if  e>ie  are  sufficiently  smooth  on  .\  then  the 

arbitrariness  of  <p  implies  (after  integration  by  parts)  that  Q(*)  *  0  in 
and  that  the  zero-Neumann  boundary  condition  on  and  the  interface  conditions 

(2.8)(c),(d)  are  fulfilled  ( ( 2 . 8) (c ) , (d )  are  the“natural"  interface  conditions  for  the 
problem  Q(#)*0) . 


A  simple  calculation  shows  that  0 (*«.)  <  0 ,  *  1 3  n  3#  *  *+  and  Q(*)  >0, 

Therefore  we  get  for  every  weak  solution  of  Q*>)  «  jj 


*! 


3ncu  3»c 


Q(# )  -  Q  (*  ^.)  *  X2div(e»7(*-#  +  )  )  -  f(x.y)(*-*  +  )  *  0 


and 

0  ( *  )  -Q(*_)  *  X  2d  i  v  ( e»7(  *-»  _  ) )  -  g  ( x  ,y  )(*-*. )  i  0 

i  *1  **i  1*1  **i 

where  f(x,y)  *  j(e  u*  *e  v*)jO,  g(x,y)  *  J  (e  u**e  v*)«0 

and  njfx.y)  U2(x*y))  is  between  ^ ( x . y)  and  As  for  the  continuity 

equations  the  maximum  principle  yields  (2. 10)(c)  . 


The  existence  statement  of  Theorem  2.1  follows  by  using  Schauder's  fixed  point  theorem 
as  in  Sank,  Jerome  and  Rose  (  1982). a 

We  remark  that  and  holds  if  a  is  empty  and  that  the  estimates  (2.10)  imply 

that  the  solution  is  unique  if  U^O  for  k»0,...,r.  Then  uaval  in  q  which  implies 
Jn«0p«0.  The  whole  device  is  in  thermal  equilibrium. 

Mock  (1974)  showed  (under  slightly  more  stringent  smoothness  assumptions  on  the 
boundary  conditions)  that  the  solution  is  unique  if  lU^I  are  sufficiently  small. 
Uniqueness  for  arbitrary  U^  cannot  be  expected  since  there  are  wellknown  devices 
(like  thyr i s to rs .see  Sze  (1981))  which  exhibit  multiple  solutions. 


We  now  show  (for  devices  without  oxide-regions)  that  any  solution  of  the  semi conduc tor 
problem  is  classical  if  the  doping  profile  D  is  piecewise  Hoderconti nuous . 

We  denote  by  C R ( ? n )  the  set  ofcritical  points  of  3q,  that  is  the  set  of  all  points 
P  €  33  for  which  either  P  €  3iij.nT?.<s  (Ohmic  contacts  and  insulating  segments  meet), 
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P  £  J  n  an. (Ohmic  contacts  and  the  oxide  boundary  meet)  or  in  OB  (P)  (where  8  (P)  is 

V  J  _  »  r 

the  open  ball  withradius  p  centred  at  P)  does  not  have  a  C  ‘  parametri sa tion  for 
any  p  >  0 . 

Theorem  2.2:  Let  1  and  assume  that  afl  is  Lipschi tzcontlnuous  and  piecewise  C* 
and  that  are  C°*  curves  for  1*1, ..,N.  Also  assume  that  0  is  defined 
on  a,  that  0£L."(o)  and  that  0|nfec°'a(n1)  for  i»0 . N  and  some 

a>0 . 

Moreover  0  is  L i pschi tzcont i nuous  in  a  neighbourhood  of  30c.  Then  every 
weak  solution  (e.u.v)  of  the  semiconductor  problem  (2.4),  (2.7)  for  which 
*  €  L*  (a  )  holds,  fulfills 

(2. 12)  (a)  (e.u.v)  £  (Cl(a  -  CS(aa)))3 

(2.12) (b)  *l3  €  C2(n1) .  i *0 . N 

and 

(2. 12)  (c)  (u.v)  €  (C2(a))2. 

Also  all  second  derivatives  of  *,u,v  are  square  integrable  over  any  subregion  of  a 
which  has  positive  distance  from  the  critical  points  of  . 

Proof:  For  a  given  weak  solution  (o*,u*,v»)  we  set 

F(x.y)  •  i2(e**(x'y)u«(x.y)  -  e"**( * *y 1 v»( x  ,y ) )  -  0(x,y). 

Since  *’,,u*,v*  €  l"(n)  we  have  F€l“(a).  **  i  s  the  unique  solution  of 

(2.13)  X2o*  .  F(x ,y) ,  (x.y)  £IJ 

subject  to  the  mixed  01 ri chi et-Neumann  conditions  (2.7)(a),  ( 2 . 7 ) (d ) ,  ( e) .  From 
Kawohl  (1980)  we  obtain  that  **€H2(a‘)  where  a'  is  any  subdomain  of  a  with  positive 
distance  from  thecritical  points  of  0  and  that  i|>  fulfills  (2.13)  and  the  boundary 
conditions  almost  everywhere  (with  respect  to  the  two  and  one  dimensional  Lebesgue 
measures  resp.).  Therefore  £ C(fl-CR( ?n) ) .  Theorem  1S.1  in  Ladyzenskaja  and 
Ural  'tseva  (  1968,  pp.203  )  implies  *•  €  C 1  ,a(  5  ’ )  and  *•  €  Cl(.T-CR(3n))  . 

The  same  theorem  implies  that  the  solutions  u*,v*  of 

d  i  v  ( e*1*^ x  ,y  ^7u  )  «  0 
and 

div(e‘**(x,y)7v)  • 0 

fulfill  u»,v»  £  C 1 ’ a ( 5  * )  and  therefore  u • , v * e C * ( n-CR( an) ) .  Theorem  6.24  in  Silbarg 
and  Trudinger  (  1  977  ,  page  6.24)  yields  u*  ,  v  *  £  C2  (  r. )  . 
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Let  n*  be  a  subdomain  of  n.  with  L 1  pschi tzcont i nuou s  boundary  and  with  a  positive 
distance  from  the  cri  tica)  points  of  .  Then  »*  ,u *  , v*  ,D*€C^  ’ a ( n" )  and  therefore 
F  €  C0’a(iT).  Since  *•  solves 

X2d*  ■  F ( x  ,y ) ,  (x.y)  6  fl“ 


we  get  »*€C2,a(n")  which  implies  **£C2(nf).  a 

*,u  and  v  may  have  s i ngu 1 ar i t i es  at  the  critical  points  of  an  and  second  derivatives 
of  w  are  discontinuous  accross  the  junction  r.  if  tO]r.  *0.  The  extension  of 
Theorem  2.2  to  MOS  devices  is  stra ight-forward  using  the  methods  of  Ladyzenskaja  and 
Ural ' tseva( 1968 ,  Chapter  3,  Section  16). 

We  remark  that  ( 2 . 10 ) ( a ) , ( b )  yield  an  a-priori  estimate  on  the  number  of  active 
carri er-pa i rs .  (2.9)(c),(d)  gives 

4 

(2.14)  n  p  ■  4  u  v  in  .1 


and  ( 2 . 10 ) ( a )  ,  ( b )  imply: 

.  .  a  /  1  VL,  J 


.15)  44exp  Vw*x')  *  np  s  a4exp  ITAiL.^  i  „ 


I  *  max  IU,  -  U  .1 
max  j  ^  k  i 

^  ma  x  I  is  the  '  argest (  i  n-absol ute  value)  voltage  applied  to  two  Ohmic  contacts. 

This  estimate  was  anticipated  by  Oe  Mari  (1968)  and  it  was  proven  by  Markowich, 
Ringhofer,  Sel berherr and  Langer  (1982  a,b)  for  the  one-dimensional  case. 

Also  estimates  for  n  and  p  follow  from  (2.9)(c).(d)  and  (2.1).  We  get  for  non-MOS 
dev  ices 


( 2 . 16 ) ( a  ) 


0.  */0. *44  exp(-lVmaxl/UT)  ( 

- - -  exp  V*  )  4  n  1 

,  ('Ly.'..)  5 

2  \  UT  / 


(2.16)15) 


.1  ' V max  1  \ 

23  eXP( - ITy  j  , 

0,  *  '/0;>4,54exp(!Vmdx!'UT) 


■>  .9  !  ^  max} 

2  5  exp  — ,v— 


D _  *•  /0_*4i  exp(  -  IVmjx  I  /Uy  ) 


I 

k  . 


r 

\  ■ 

fc- 

F 


If  0_<0,  0+ >  0  then  the  lower  bounds  in  (2.16)  ere  0(6*)  as  6  -*0,  the  upper 
bounds  are  0(1)  (for  fl xed  lvmax0 • 

There  fs  numerical  evidence  (see  Markowlch,  Rlnghofer,  Sel berherr  and  Langer 
(1982  a.b))  that  the  estimates  (2.16)  are  not  sharp.  The  factors 

exp  ( — •  *xp( — tj*  c*n  probably  be  omitted. 
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3.  ASYMPTOTIC  EXPANSIONS 


We  are  now  concerned  with  the  asymptotic  behavior  of  the  solutions  of  (2.4),(2.S), 
(2.7) ,  (2.8)  as  x  -0+. 

When  we  set  x»0  In  (2.4)(a)  we  see  that  * ( x • 0 )  has  to  be  discontinuous  at  pn,  nn 
and  pp  junctions  ( s 1 nceO  1$ d 1 scontl nuous  at  these  junctions)  assuming  that  u(x«0) 
and  v(x»0)  are  continuous  In  n  (note  that  (2.4)(b),(c)  only  admit  weak  solution  In 
C(n)  If  ii  E  L"(0)).  Therefore,  standard  singular  perturbation  theory  Implies  that  we 
have  to  expect  an  internal  layer  in  «  (that  is  a  region  of  fast  variation  of  *  )  at 
these  junctions.  Also  boundary  layers  occur  at  an  if  the  reduced  solutions  (i.e.  the 
solutions  of  (2.4)  with  x«0)  do  not  fulfill  the  boundary  conditions. 

for  the  following  analysis  we  assume  that  the  profile  0  is  discontinuous  accross  only 
one  (open)  C“-curve  r  which  splits  n  into  two  connected  subdomains  n  +  and  n_.  Also 
we  assume  that  0!-  £C®,a(n  +  ),  €  C^'“(n_)  for  some  a>0  and  [D)r  *0.  Also 

?  n  an.  ■  {  >.  0  does  not  change  sign  in  n+  and  in  n_. 

We  denote  by  t(x,y)  the  or i ented d i stance  of  (x,y)  from  r,  that  means 

t>0  in  n+  and  t<0  in  n_ .  s(x,y)  •(s.(x,y),s.(x,y))  is  the  point  on  r  which  is  closest 
to  ( x ,  y )  (s  is  unique  in  a  sufficiently  small  strip  about  r).  Similarly  r(x,y)  >0 
denotes  the  distance  of  (x,y)  to  3n  (an  fulfills  the  assumptions  of  Theorem  2.1)  and 
<4(x«y)”(Q1(*.y)  •  q2(  x  ,y )  )  denotes  the  point  on  an  closest  to  (x,y).  Note  that 
7r'an’"*  and  7t I r  ’S  the  unit-normal  vector  of  r  pointing  into  n+. 

For  a  function  f  defined  on  a  (or  n)  we  set 

fr(t.S)  ■  f(x,y) 


and 


fa(r,q)  ■  f ( x , y ) 

in  neighbourhoods  of  r  and  an  where  s  and  q  resp.  are  unique. 


'V e  define  for  some  s€T 


f  r ( Or ,  s )  :  *  lim  f(a,b) 
(a ,b)  -  s 
( a  ,  b )  €  n  + 

f ‘ ( 0-  ,  s )  :  *  lim  f(a,b) 
(a , b )— s 
(a ,b)en_ 
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(assuming  that  the  limits  exist). 


We  also  set 


'3S1 

3*2*1 

( a<3> 

3q?'i 

3X 

3  X 

.  sT^Ty)  * 

TT 

TT 

3S1 

V  3y 

3 s2  j 

ay  J 

3qj 

3q2 

ayj 

f3f 

(if*} 

3  ft 

*  V*  * 

To7| 
11* ! 

and  remark  that 


•  3 


holds  (the  superscript  T  denotes  transposl tion) . 

Following  standard  singular  perturbation  theory  we  make  the  ansatz': 


/f(x.y.A)' 

(x.y)' 

5,-(*.y) 

.^(x.y); 

'$.(11^1,  s(x.y))' 

q(x  ,y ) )  | 

u(x.y.A) 

U(x,y,A)J 

~  I  A1 
1-0 

♦ 

s(x  ,y) ) 
s(  x  ,y) )  j 

* 

q  (  x  ,y ) ) 
7i(li*xJd,  q(x.y))^ 

where  the  functions  marked  with  are  independent  of  A,  the  functions  marked  with 
'a'  are  defined  on  (-<»,*)  *  r  and  decay  to  zero  as  r-  *«•(  i  nternal  layer  terms). the 
functions  marked  with  are  defined  on  (0,®)  x  a  12  and  decay  to  zero  )))■[-■ 
boundary  layer  terms). 


We  insert  (3.1)  Into  (2-4),  (2-5),  (2.7),  (2.8)  and  obtain  equations  for  the  i-th  teunln 
the  series  (3.1)  by  comparing  coefficients  of  a  .  We  start  with 


A)  The  Reduced  Problem 

Evaluation  of  (2.4)  away  from  r  and  Ail  and  comparing  0(1)  terms  gives  (after  dropping 
the  index  0)  the  zeroth  order  reduced  problem 


( 3 . 2) (a) 

0  .  a2e'*'u  -  i2e"*v 

( 3. 2) ( b) 

d i v ( e*cC)-  0 

( 3 .  2 )  ( c ) 

div(e"vov) *  o 

0(x,y)  ) 

t 


1  (x  ,y)  €n  -  r 

i 

i 


In  the  context  of  semiconductor  device  physics  this  problem  is  referred  to  as 
'zero-space  charge  approximation'. 
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8y  Investigating  the  Interna!  and  boundary  layer  problems  we  will  supplement  (3.2) 
by  Interface  and  boundary  conditions. 


8 )  The  Internal  Layer  Problem 

We  evaluate  (2.4) ( b )  close  to  r  but  away  from  so  and  compare  0 ( x ” ^ )  terms  (after 
carrying  out  the  differentiations) . 


Thi s  yields: 


u„  + 

Ott 


U  tl 
Ot'O  T 


0 


( subscripts^*  denote  differentiation  with  respect  to  t)  and  therefore 

“0T  •  *(s)e  0T  ■  F  rom  uQ  (  t  ,  s )  » 0  as  t-*»we  get  u0  ■  0 .  Similarly  vQ  ■  0  follows,  u  and 

v  do  not  have  zeroth  order  internal  layers. 

A  A  A  A  A  a 

For  the  following  we  set  u  ■  u  ^ ,  v  *  v  ^  and  ^  ^  . 

Comparing  0(X  )  terms  in  (2.4){b)  (again  after  differentiation)  and  evaluating 

close  to  r  but  in  n+  gives 


2rT  ♦  *T(Gtr(°+,s)  +  uf)  •  °.  r>0 

(subscripts  t  denote  differentiation  with  respect  to  t)  and  evaluation  close  to  r  but 

inn: 


Integration  yields 


0  ,  t  <  0 . 


(3.3) 


2t(t.s) 


3^(0+.s)(e**(t,s)-l) ,  t  >  o 
r  l 

^  G^(0-.s)(e'®(T,s)-l) ,  t  <  0 


Proceeding  analogously  with  ( 2 . 4 )  ( c )  gives 

i[(0+,s)(e*(T,s)-l),  x  >  0 
^(0-,s)(e^(t,s)-l)  ,  r  <  0 


(3.4)  vt(t.s)  *  { 


We  used 


that  $ , u  , v 


and  their  t -der i va ti ves  vanish  at  ***<■>  for  all  s€r. 


The  internal-layer  problem  for  I  is  obtained  by  evaluating  (2.4) ( a )  close  to  r  (but 
away  from  )n)  and  by  comparing  0(1)  coefficients: 
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(3 . 5) (a) 


l42eir(0+.s)4u-r{0>s)  .52#-*r(0*.s)-«v-r(0.fS) 


-  0  (0+,s) , 


t  >  0 


(3. S) (b) 


•  „  .  42#il(0-.s^u-[o.,j).42#-;r(0-,s,^;r{0.>s) 

-  0r(0-,s),  t  <0 


Interface  condi  tion  for  ( 3 . 5)  (a ),(b)  are  derived  by  using  that  *  €  C^n)  (see  Theorem 
2.2),  which  Implies  *r ( 0+,s ) «*r (0- ,s )  and  *» (0+,s)«*[{0- ,s ) .  Inserting  (3.1)  into 
these  relations  and  comparing  0(1)  and  Q(x*^)  coefficients  yields 


( 3 . 5 )  (c )  *(0+,s)  - C(O-.s)  «  ir( 0-  ,  s )  -  ir(0+,s) 

(3. 5)  (d)  JT(0+,s)  •  JT(0-,s). 

( 3 .  5 )  ( e )  *(+<», s)  ■  i(— ,s)  *  0 

A 

(3.3)  -  (3.5)  are  supposed  to  hold  for  all  s€r.  *  is  discontinuous  at  t*0  for  all 
s  €  r  since  Ml.  •  0. 


Theorem  2.2  also  provides  interface  condi ti ons  Tor  the  reduced  problem  (3.2)  since  it 
implies  that  ur(0*,s)  «  u‘(0-,s),  vr(0+,s)  *  vr(0-,s)  and  u^(0+,s)  *  u£(0-,s), 
v£(0+,s)  *  Vj(O-.s).  Inserting  the  expansions 

(  3 . 6 )  (  a  )  u(x.y.X) -u(x.y)  ♦  xu(£,s)  +  ... 

( 3 ■ 6 ) ( b )  v (x ,y , A )  ~  v(x ,y )  +  a$(£,s)  ♦  ... 

(the  dots  denotea  power  series  in  \  stating  with  the  0 ( x *• )  term)  into  these 
relations  and  comparing  0(1)  terms  yields 


( 3 . 7 ) ( a )  (G)r  *  0,  C v ] r  •  0 


and  using  (3.3),  (3.4) 

G^(0+.s)e'^{0*,s) 

;[(0>.s)eJ(0+-s> 


u£{0-,s)e'*^°"’s^ 

5[(0-,s)e*<°-'s) 


(3.5)(c)  implies 


r  Tr 

u>  ( 0  + , s ) e ^ 

r  7r 

v^(0-*-,s)e_t' 


( 0+ ,  s ) 
( 0* ,  s ) 


(0-,s)eir(0‘,s> 

(0-,s)e‘*r(0*,s). 
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So  we  obtain  the  interface  conditions  for  the  reduced  problem 


( 3 . 7 )  ( b)  [e*7u'7tlr  *  0,  [e**7;.7t)r  •  0  . 

since  •  7f  (x  ,y )  ■  7t(x  ,y ) !  ( x  y }  m^T  ho1  ds . 

C)  Ohmic  Contacts 

A  straight'f orward  calculation  shows  that  the  Oirichlet  boundary  conditions  on  3Qj. 
for  i|« ,  u  and  v  fulfill  (3.2)(a).  Therefore  we  do  not  expect  boundary  layers  at  3(ic- 
T  hat  means 

(3.8)  *0  (<>.<?)■  7(3  (0.d)"~o(0.q)»O  for  b  >  0  and  qUQj.. 

To  proove  (3.8)  one  has  to  proceed  as  Markowich,  Ringhofer,  Selberherr  and  Langer 
(  1982  a,b)  did  for  the  one  dimensional  semiconductor  problem. 


(3.8)  implies  that  we  have  to  impose  the  same  Oirichlet  boundary  conditions  for  the 
reduced  problem  (3.2)  as  for  the  full  singularly  perturbed  problem: 


(3.9) 

where  ti,u,v 


* 1  3ac  *  * ' 3nc  * 

on  sri.  are  given 


by  (2.7). 


anr 


1  an- 


0)  Insulating  Segments 

We  assume  that  0  is  differentiable  in  a  neighbourhood  of  3(1  and  that 

(3.10)  70-7 :  .o 

i  s 

holds.  Differentiating  (3.2)(a)  gives 

(3.11)  0  *  S2(e*u+e'^v)7j  +  52(e*7g-e'll'7v)  -  70 
Since  e^u+e  *v  >  0  in  ii  we  get 

(3.12) (a)  7j-7‘  .  o 

1  s 

i  f 

(3.12) (b)  7u'n;?i7  •  7v.ni3  a  •  0 


holds.  (3.2)(a)  is  compatible  with  the  zero  Neumann  conditions  for  *,u,v  on  jfl  ^  $  and 
we  get 


(3.13) 


*0(o,q)  ■  u0(o,q)  ■  v0(o,q)  ■  0  for  p>  0,  q  €  312.  $ 


as  well  as 


(3.14)  (0  ,q )  ■  Uj(e  ,q)  ■  7j  (e  ,q)  ■  0  for  p>0,  q  €  siij  s  . 

No  zeroth  and  fi  rst  order  layers  occur  at  312<S. 

If  we  did  not  assume  (3.10)  then  (3.14)  and  ( 3 . 12) ( a)  would  not  hold  (of  course 
(3.13)  would  still  be  valid). 

( 3 . 12 ) ( a ) , ( b)  define  homogenous  Neumann  boundary  conditions  for  (3.2)  on  3fl^  . 

E)  Oxide-Semiconductor  Interface 

As  for  junctions  we  get  ufl  ■  7fl  ■  0  and  setting  u»up  7  ■  7j ,  7  »  we  obtain 
(3.  IS) (a)  ua(o.q)  ■  G ’ ( 0 , q ) ( e ‘ < 0  * q > - 1 )  .  o  >  0  .  q  €  3nQS 

(3.15) (b)  vs(o.q)  •  ;®(0,q)(e,l(‘1,q).l).  «>0.  q£»JS' 

Since  we  set 

( 3 .  16) ( a )  u  (  x ,y  ,  X )  ~u(x,y)  *  Xu(£,s)  +  Xu(£,q)  *  ... 

(3.  16) (b)  v ( x  ,y  ,  X )  ~  v ( x  ,y)  +  Xv(y,s)  ♦  x7(£,q)  +  ••• 

1 

(where  the  dots  denote  a  power  series  in  x  starting  with  the  0 ( x  )  term),  we  get 
by  inserting  (3.16)  into  (2.8)(e)  and  by  comparing  0(1)  coefficients 


0  -  ( Vu  -  n )  ( q ) 

♦  a*(0,q)  (e”l'(0*q)-l)  . 

q  €  312 

0  • ( V» • n) ( q ) 

♦  v*(0.q)(e^0><,)-l). 

q  e  312 

Therefore 

0  «  u’(0,q)  «  v’(0,q)  ,  q  £  3S2QS 

and  7(3, q)  ■  7( o  ,q )  ■  0  for  s  >0,  q  £  3120J .  Boundary  conditions  for  u,v  on  5nos  are 


(3.17) 

If 


7u  • 


312 


os 


7v  ’ 


3W 


OS 


0. 


(3.18) 


70.  ni  *  0 

9%s 
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I 


holds  then  (3.17)  Implies  (as  for  insulating  segments) 

(3.19)  7*-n|  .  0. 

3nos 

By  evaluating  close  to  JflQS  but  away  from  r  we  obtain  the  boundary  layer  equation  for 

«|  which  is  analogous  to  the  interface  problem  (3.5): 

*nos 

( 3 . 20 )  ( a )  *BD  -  42ei'3(0’<,)+*u3(0,q)  -  6  V** !  0  *q)  '*J3(0  ,q)  -  03(0.q).  p>0.  q  €  jqqs 

(3.20)  (b)  *(-.q)  *0,  q  £  300$. 

To obtai  n the  boundary  condition  for  (3.20)  at  p*0  we  solve  Laplace's  equation  in  the 
oxide. 

Let  G(x,y,C,n)  denote  the  Green’s  function  (seeProtter  and  Weinberger  (1967),  Chapter 
2,  Section  7)  of  the  problem 


o<p  «  f  i  n  t 


<fli 


3e, 


g ,  7-p .  c ; 


«  h,  7(f>-  cl  .  k 

His  ■  3nos 


(note  that  e'  ) ■ 

3nos  3nos 


Then,  since  *»*G+«p  fulfills 


Oip  «  0  in 


3a, 


0,  7<e-;| 


3  # 


is 


OS  0 


3(5 


OS 


we  get 


J  G( x ,y  ,t ,n)7*(t ,n) • «(S ,n)d(C ,n) 

0  J0os 

This  gives  the  boundary  condition  (for  the  singularly  perturbed  problem  in  (5): 


(3.21 


1  3(5 


OS 


^  5^0SG<X’y’"’n) 1 (x.y)€anos 


7*(? ,n) • n(? ,n)d(5 ,n)  *  *G- 


For  the  following  we  assume  that  if  r  hits  3^os  (which  in  fact  happens  in  MOS-techno- 
logy)  then  F  is  perpend i cu 1  a r  to  Tnos,  that  means 

(3.22)  7t(S)’7r(S)  *  0  for  (SI  ” f  n  s . 

We  Insert  the  expansion 

(3.23)  *(x,y  ,  \  )  ~  I  (  x  ,  y )  ♦  *(i,S)  *•  *(^,q)  +  ... 
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I 


(the  dots  denote  a  power  series  starting  with  the  0(x)  term  whose  coefficients  have 
the  same  form  as  the  0(1)  terms)  into  (3.21).  Assuming  that 


(3.24)  Ct  <  I *G I x  <  C2 

where  C^.Cj  are  independent  of  x,  we  obtain  by  comparing  0(X**)  terms  and  by  using 
(3.22): 

(3.25)  J  S(x,yfC.n)|(x  yj€Jfl^^-*o(0.q(e,n))d{5,n)  •  X^Yg- 

This  integral  equation  is  uniquely  soluble  (for*  ( 0 ,  q ) )  since  *  (O.qH'fw-;)  (q) 

P  P 

where  w  is  the  (unique)  solution  of 


w|3«,  * 


■  0, 


i  s 


1  3P. 


Obviously  w  €  c"(i-CR(ae))  holds. 


If  >  is  a  rectangle  as  in  Figure  1  (which  is  a  common  assumption  in  MOS-model 1 ing) 
with  d«dist  (30os»a*c)  then  w  is  a  *iftear  function  and^w"t|,-  *  -x-^*..  In  this 
case  os  Es 

(3.26)  »0(Q.q)  *  -X'jjr"  *2  • 

No  zeroth  order  layer  occurs  at  the  oxide-semiconductor  interface  if  the  right  hand 
side  of  (3.26)  is  not  0(1)  as  \-0+. 

(3.25)  (or  equivalently  (3.26)  if  *  is  a  rectangle)  provides  the  missing  boundary 
condition  for  the  interface  layer  problem  (3.20). 

Equations  for  the  higher  order  terms  of  the  expansion  (.’.1)  *,1  be  derived  in  a 

analogous  way. 


When  the  asymptotics  of  *,u,v  are  known  then  expansions  of  n, p  and  of  Jn,  Jp  as  given 
by  (2.9)  can  easily  be  derived,  We  get  from  (3.1) 

(3.27)  n(x,y,x)  »  n ( * , y )  ♦  fl(j.s)  +  n(j,q)  +  ... 

(3.28)  (b)  D( x ,y , x )  ■  p ( x ,y )  ♦  S(j,s)  +  P(£.q)  ♦  ••• 

(3.16)  and  (3.23)  imply 


(3.29)(a) 


2  ii- 

«  Vu, 


n(t,s) 


«2e*  CO+.s) (eqi(T ,s)_1  )Qr(0,s) ,  t  >  0 
42e+r(o-,s)(e*(T,s)_1)ir(0>s)>  T  ,  0 
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( 3 . 29 )  ( b )  ri  ( 0  ,  q  )  ■  0  for  a  >  0 ,  q  £  i(l|.  U  5Bj 

(3 . 29)  (c)  n(a,q)  »  i26t'5(°lt’)(e"::(p,q)-l),  o  >0.  q  €  3nos 

(3.29)  (d)  p  .  jZe**5 

.  .  A  — 

and  similar  expressions  for  p  and  p. 


Differentiating  (3.16)  and  using  (3.3),  (3.4)  gives 


( 3 . 30 ) (a )  Jn(x,y,X)  .  J„(x.y)  *  3n(|.*)  *  Jn(C,q)  ♦  ... 


( 3 . 30) ( b) 

V 

(x.y.A)  •  3p ( x ,y )  ♦  Jp(j.s)  ♦  Op(j.q)  +  ... 

with 

( 3 . 31 ) ( a) 

^n 

*  s2e*wa,  (b)  Jp  —«V*7» 

(3,32) (a) 

A 

Jn< 

'i2eir(0+,s)(e:(r,s).1)  isl^l,  .7  jr(0r,s),  r>0 

i  t  .  s )  »  \  "  M*,y)!(x.y).s  5 

!.  3 2  e  J  ‘  (0-.s)(ef(T,s)-i)  istXjjM  .7G:(0-,S),  t<0 

1  ,y,!(x,y)*s  5 

( 3 • 32) ( b) 

Jn< 

!  o  ,q  )  ■  0  for  o  >  0  ,  q  €  30c  U  3(1,-  s 

( 3 ,32 ) ( c  ) 

V 

:  o  ,q )  *  52e*i(0’q)(e*(o’q)-l)^|lilj  j  ^  t7qG3(0.q) .  q£;qJS 

A  ~ 

Analogous  expressions  hold  for  Jp  and  Jp. 

Since  (T[^T7T)T'7tir  *  (7nbyy)T7rl3r.  "  0  holrfs  get 


(3.33! 


■7tl 


7r  i 


3(1 


■  w 


OS 


Jp-Vr| 


3(1 


■  0 


OS 


Therefore  the  current  density  components  perpend i cu lar to  the  j unc ti ons ( to  the  semi¬ 
conductor  oxide  interface)  have  no  zeroth  order  layers  while  the  current  density 
components  parallel  to  the  junctions (to  the  semiconductor-oxide  interface)  may  very 
well  have  zeroth  order  layers.  This  phenomen  is  illustrated  by  the  MOS-transi stor 
simulation  performed  by  Selberherr  ( 1980) - 


So  'ar  we  only  considered  one  curve  of  discontinuity  of  0.  Generally,  an  internal 

layer  in  p  occurs  at  each  junction  and  each  layer-term  * .  fulfills 

the  corresponding  layer  prob’em  (3.5)  (with  rv-local  coordinates  (s^,t^)). 


4.  Existence  of  the  reduced  solution  and  of  the  layer  solutions 


A.  Reduced  Solution 

We  now  assume  that  we  have  N  C* -junctions  which  are  as  in  Section  1.  [D]ri  *0 

holds  for  i»l,...,N.  We  denote  the  normalvector  to  by 
Then  the  interface  conditions  (3.7)(a),(b)  read 

( 4 . 1 ) ( a )  [u]  *0,  [ v ]  -0,  i *  I . N 

"  i  ri 

(  4  . 1 )  (  b)  [  e*7u  '  tj  ]  r  *0.  [e"*7v-  t.  J  «Q,  i«l . N. 

1  ri  1  ri 

Moreover  an  fulfills  the  assumptions  of  Theorem  2.2.  We  prove: 

Theorem  4.1  Assume  that 


0  €  C(ni),  i-1 , 


holds  and  that  0,  r^,  an  fulfill  Che  assumptions  given  above.  Then  the 
reduced  problem  (3. 2), (3. 9),  ( 3  .  1 2  )  ( b )  ,  ( 3 . 1 7 )  ,  ( 4 . 1 )  has  a  weak  solution 
(u.v)  €  (Hl(n)  n  L“(n) )2,  JeL“(n).  Every  weak  solution  fulfills  the 
estimates  (2. 10) (a) »(b)  and 

(4.3)  *  i(x.y)  s  (x.y)  e  n 

Mere  7. ,  are  defined  in  ( 2  . 11 )  ( a ) ,  (  b ) ) .  Also  every  weak  solution 
satisfies 

( 4 . 4 )  (a  )  l  €  C ( 5 i  -  CR( an) ) ,  (b)  ( G ,J )  e  i C ( 5  -  CR( an ) ) ) 2  fori-0 . N. 

Proof ;  A  weak  solution  of  the  reduced  problem  is  given  by  a  triple  (i,u,v)with  (  u  ,v  )fi(  H  (il) )  ^ 
J  E  L°"(n)  which  solves  ( 3 . 2 ) ( a )  pointwise  (almost  everywhere),  which  assures 
the  Oirichlet  boundary  conditions  on  an^  in  the  sense  of  (H*(n))2  and  which 
fulfills 

(4.5) (a)  Lj(*,u,<p)  ■  f  e^vu '  Pujdx^y  *  0 

n 

(4.5) (b)  L  2  ( •  v  ,tp)  ■  Je'^vv  •  vgdxdy  «  0  for  all  $  e  h*(  n  u  an  ^  u  anos ) . 

n 

(If  <i  ,u  ,v  and  are  sufficiently  smooth,  then  (4.5)(a),(b)  yields  (  3 . 2  )  {  b  )  , 

(c),  (3.12)(b),  (3.17)  and  (4.1)  using  integration  by  parts.) 

The  maximum  principle  for  weak  solutionsof  linear  equations  in  divergence 
form  immediately  implies  that  u.v  fulfill  (2.1)  (a),(b)  if 
(3.2)(a)  gives 


I 


(4-6) 


J  (  u  ,  v  )  *  In 


L  Z4‘u  J 

and  (4.3)  Is  immediate.  (4.4) ( b )  follows  from  a  posteriori  regularity 
theorem  (see  lady zenskaja  and  Ural  1 tseva(  1968) ,  Chapter  3,  Section  14)  and 
(4.4) ( a )  Is  implied  by  (4.6). 

To  establish  the  existence  statement  of  the  theorem  we  employ  Schauder's 
fixed  point  theorem  similarly  to  the  existence  proof  of  Bank,  Jerome 
and  Rose  (1982)  for  the  singularly  perturbed  problem  (see  Theorem  2.1).  The 
existence  proof  for  the  reduced  problem  isslightly  more  difficult  because 
the  regularising  property  of  Poisson's  equation  is  lost. 

We  define  the  following  set  M: 

(4.7)  M«((u,v)£( l“(n ))2iu_*u<uv,  v_Sv«v+) 

2  2  2  2 

and  regard  M  as  closed  and  convex  set  in  (L  (a))  .  A  mapping  T:  M-»(L  n)) 
is  def i ned  as  follows: 


(4.8) 

T(Vvo>  *  <ul*vl> 

where  Uj.Vj  are  the  weak  solution  of 

( 4 . 9) ( a  ) 

div(e*(u°,vo)’7u1)  *  0 

( 4 . 9 )  (  b ) 

div(e*^u°,vo^7v1)  -  0 

subject 

to  the  already  specified  boundary  conditions.  T  is  well  defined  since  (4.9) 
(a),(b)  has  a  unique  weak  solution  ( u j  ,v ^ )  6  ( (fl ) )  .  Clearly  every  fixed 
point  of  T  s  a  weak  solution  of  the  reduced  problem. 

The  a  priori  estimates  ( 2 . 10) ( a )  ,  ( b)  imply  that  T : M  ■*  M(  since  they  are  inde¬ 
pendent  of  the  mul t i pi i ers of  vU)7v  in  (4.9)).  The  wel 1 posedness  of  (4.9)  in 
(  H * ( )f i nduces  the  boundedness  of  Range(T)  in  (H*(3))^,  therefore  Range(T) 
is  relatively  compact  in  (  L  2  ( .1 )  )  ^  . 

To  show  that  T  6C(M)  we  take  a  sequence  (u  ,v  )  €M,  such  that 
1 2  ( " ) 

u  j  *  n  -  ® 

n 

i2(.u 

v  -•  v  • ,  n-*«» 

n 

holds.  We  will  show  that  every  subseouence  of  T(j  #vn)  has  a  subsequence  which 
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converges  to  T(u*,v*) 
(L2(n) )2) . 


in  ( L 2 ( n ) ) 2  (which  implies  T(un,vn) 


T ( u * , v* )  in 


Since  the  range  of  T  is  relatively  compact  in  M  and  since  M  is  closed  we 
conclude  that  every  subsequence  of  T(un,vn)  has  a  convergent  subsequence  (with 
limit  in  M ) .  let  T ( u n  ,vn  )  be  such  a  convergent  subsequence  (of  a  subsequence 
of  T ( u  , v  ) )  and  let  k  k  T(u„  ,v„  )  -  (a.b)  €  M  in  ( L2 ( n ) ) 2  as  k 

n  n  n  ^  n . 

Obviously  l“'  * 

L  2  ( )  _ 

*(un  ,vn  )  ■»  *(u*,v»)  as  k-«» 

holds.  Therefore  there  is  a  subsequence  i ( u n  ,vn  )  such  that 

kj  kj 

*(V  >V  > 

e  J  J  ,  ei(u*.v») 

j  -- 

•*(%  •  vn  ) 

e  #-i(u*,v«) 

almost  everywhere  pointwise  in  a  . 

Theorem  5.4  in  Ladyzenskaja  and  Ural  1 tseva ( 1968,  Chapter  3,  Section  5)  implies 
that  the  weak  solutionsof 

-V  > 

k  .  k  . 

div(e  J  J  tru )  »  0 

";(uV,vV! 

d  i  v  ( e  'l  ^7v)*0 

subject  to  the  boundary  and  interface  conditions  converge  in  H 1  ( n )  as  j 
to  the  solutionsof  the  limiting  problems 

div  (e*(u  *'v  *^7u)  »  div  (  g  u  *  ■ v  * )  ?v  )  *  0 

(subject  to  the  boundary  and  interface  conditions). 

This  yields  T(u„  ,vn  )  -  T( u *, v «)  i n  ( L 2 ( a ))2  as  j  ■» «•  and  (a.b)  »  T(u*,v»). 
n  k  n  k 

Therefore  T  is  j  j  a cont i nuous  and  compact mappi ng  from  M  into  M  and 
Schauder's  fixed  point  Theorem  implies  the  existence  of  a  fixed  point  of  T 
in  H.  a 


Note  that  I  has  jumpd i scont i nu i t i es  along  r ^  ,  i*l,...,N. 

For  the  carrier  densities  n  and  p  we  get  from  ( 3 . 29 ) ( a )  ,  ( 4 . 6 ) 
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n  * 


0  +  '/o2»4i4uv 


P  * 


s  *  { 


0(4 


I  v. 


expt — 07" 


■0  +  sV 

+ 

X* 

o. 

C  1 
<  1 

2 

;  on  u , 

v  yield 

-) I1  n  ai 

^  01  = 

i 

in 

if  DIS. 

in  L  (SI) 


and 


0  ( 4 4  exp('  /"'ax' 


n - })  in  j.  if  0  1  ;  >0 

UT  1  i 

A  1  V_.  J 

-0  ♦  0(44  exp(-u^-a—  ))  in  3i  if  0|;  <0 


in  L  ( a)  . 


If 


IV 


44exp( — -jX—  )  is  small  (‘low  injection  condition1)  then  the  electron  density  in  the 
zero  space^charge  a pprox i ma t i on  is  close  to  the  doping  profile  in  n-regions  and 
close  to  zero  in  p  regions  while  the  hole  density  is  close  to  zero  in  n  regions  and 
close  to  the  negative  doping  profile  in  p  regions,  n  regions  are  depleted  of  holes 
and  p-regions  are  depleted  of  electrons. 

We  now  investigate  the  behaviour  of  I.u.v  as  4  -0.  We  assume  that  there  is  exactly 

one  Ohmic  contact  C  •  in  every  .1  ^  for  i  *0 . N  and  that  no  two  n  (and  no  two  p) 

regions  have  a  joint  boundary,  that  means 


(4.10) 

sgn  0i 

'•  i  - 1 

sgn  0 ' ;  , 

holds. 

We  define 

the  functions 

e*  &T  in 

pj  if  D 

(4-U) 

=o  *  { 

“o,  in  2 

i  if  oi5, 

where  u 

0  1 

solves 

la) 

d  i  v  (  - 

* 

0  in  J. 

(i>) 

(4.12) 

7U07 

r‘i?is 

u  ,4os>  = 

(c) 

%lc. 

'Jr 

=  e 

(d) 

^o'u 

=  0 

=  0 
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I 


and 


(4.13)  •  vfl  - 


{ 


u . 

1 

St 

e  in  If  Dig  <  0 

in  if  0  I g  >  0 


01 


where  v  Qi  solves 


■U- 


a)  div  ( g Q )  -  0  in  n. 


(4.14) 


;b)  7vo',l'an1n(inisU3Qos)  *  0 


(c)  v 


o  '  C , 


^o^n,-aa*  3  ’ 


(p-region) 

(n-region) 


The  condition  (4.10)  implies  that  (4.12)(d),  (4.14)(d)  give  Dirichlet  boundary 
conditions  for  u^,  vQi  resj.  by  using  (4.12).  The  problems  (4.12),  (4.13)  have 
unique  solutions  u  vflj  €H  (a()  (see  Kawohl  ( 1978 ) ) assumi ng  that 


(4.15)  Oi-  e  C 1 ’ a ( a < ) ,  i *0 , . . .  ,N . 

“i 

We  prove 


Theorem  4.2.  let  the  above  mentioned  assumption  on  a,0  hold  and  assume  that  fulfills 
the  assumption  of  Theorem  4.1.  Also  assume  that 


(4.15) 


ex.p  ( 4max 

i 


) ( exp ( 


IVm,J 

™i2-)  -  1) 


<  V 


where  v and  5  are  sufficiently  small  (depending  on  D,  nQ  ...  .7^  only)  and  that 
7%i  e(L"(ni ))  ^  ,  7^oi  ))2h°l  ds  .  Then  the 

reduced  problem  (3.2),  (3.9),  ( 3 . 1 2  )  ( b )  ,  (3.17),  (4.1)  has  a  locally 
unique  solution  (u,v)  e(H^(q))^  which  fulfills 


(4. 17)(a) 


(«. 17) (c) 


i n  L .!■}. 


u  ■  uQ  +  0 (y ) 
and 


(b)  v  »  vo*0(y) 


ii 


In(^)  ♦ 

{,  3 

Int^)  ♦ 


The  estimates 


j^+  0  (y  ♦  s  )  i  rv 
T 

IJ  . 

n— -  +  0(v  +  5  ^ )  in 

i 

are  uniform  in 


0!c  >0 

“M 

o;-  <  o 

(as  long  as  (4.16)  holds). 


Prop  f ;  We  assume  that  there  are  only  two  differently  doped  regions  .<0,,7j  separated 
r  and  that  D'c  <0,  O'-  >0  holds.  The  generalisation  (under  the  given 


o 


assumptions  on  o  and  0)  to  more  than  two  differently  doped  regions  is  straight 
f orward . 

We  set 


u  *  uo  +  £u-  *  '  vo  f  £v 


*  *  S  +  { 


0  U1 

*(Tf  in 


t2  uo 

1  n  nc 


l 


1  n 


and  compute  £  from  (3.2)(a) 

i Q<Vo^+464uv 


(4.18) 


£*  *  E*(EutEV5) 


'1 


j20eio(  )u 
20exp  (-JS)v 


t  1  n 


i  n  n 


1 


in  .2  _ 


[  -0+iuM5uv 


Then  we  rewrite  ( 3 . 2)  ( t>) ,  (c ) ,  (4.1)(a),  (b)  as 


E, 


(4.19) 


d i v  ( Oe  %E  )  ■  0  inn, 

/  °o~ul  ,  E  \ 

div  U4e  0^~"  t^V)  e  ^Zu*7u0''J  *  0  in  *5o 
E  UoEUl 

CEuJ,  *  0.  °e  "J^Tt:r(I)  *i4e_^T"(-  C)<7Eu*7V'rt  r(0).»^Jaisu9 


su3no 


where  f  | .  (a,b)»  lim  fvx.y)  for  (a  ,b)  €  r  and 

‘(U  (*.y)-(a,b) 

(x.y)£n. 


f  VU1 

!  — 0 —  .£ 

div  , J4e  T  |t  TEy+7vo; J  •  0  in  a 


(4.20) 


-E„ 


div  (-Oe  7Ey)  *  0  i  n  n0 


(£„ 


-E 


VU1 


0,  -Oe  ;'?£v7t 


(C) 


5‘e  ^  (b)«  E*(^v*0)«:r{i)7Ew-3iJnj#UJ 


We  define  the  operator  Ps  :  (H*(  .2  u  Jfl,  s  u  3S2qs  ) )  2  -  ( ( H*(  3  J  3  0  .  s  u  3fiQS  ) ) 2 ) ' 


(the  superscript 


denotes  the  dual  space)  by 


•Vipjdxdy 


(4,21)  F^.E^.^) 


J 

al 

U  -U. 

4  “Hr —  1  ^4 ( Z1 ’  Z2 *4 ^ 

5  8  DT  I  (rjj>«  (7£1+7u0)'7«»1dxdy 

ao 


♦  S  (-O)e 


•E*(E1.E2.6) 


*E, 


VU1 


*  ,«t“BT  rl  ,‘E*(E1-E2 


•  7<pjdxdy 

.4) 

(7E2+7v0)7<Cjdxdy 


for  all  («pl  ,<o2 )  €  (H*(.l  u  3a.$  u  Jfios)  >2. 


3d, 


3fl, 


yields  F5(EU  ,EV)  *  0, 


that  means  F,(E..,EU)  *  0  is  the  weak  formulation  of  (4.19),  (4.20).  In  the  sequel 

.  a  u  v  , 

we  equip  (H‘(a  u  Jilj$  u  3il„s))  with  the  weighted  norm 

VU1 

II  (Fj  ,F2)1|  2  *  /  I7f.  i2dxdy  ♦  i4e  ^  JlvF^dxdy 

4,1  nx  a„ 

V«1 

*  /  |7F.|2dxdy  *  J4e  J  l7F?lZdxdy 

ao  nl 

and  corresponding  scalar  product.  The  dual  space  is  equipped  with  the  induced  functional 
norm  II  •  I!  *  ^.)  3ecause  Of  (4.11)  -(4.14)  and  (4.18)  F"  Q  ( 0  , 0 )  »0  holds  and 


(4.21)  IIF^O.Oll*^  *  °(<' 


where  < 
given  by 


IU  ! 

v/' exp(  3  max -g— ) .  The  bilinear  form  obtained  by  linearizing  Fq  at  (0,0)  is 
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J 


L(3l'J2,'W  ’  ,j  ‘^idxdy  ♦  54e  J  (  qj)  (?o t  •  7^ 

*4  1  ^  _ 


*  e  a27(^o  ’  )dxdy 


/  (-0)7o2-7<o2dxdy  +  44e  t  l(?o2-V<e7 

0  3  ^ U  *  * 

ui 

*  e^jjVv  ■  7(«2  )dxdy . 


To  show  that  L  is  coercive  we  estimate 


I  o,7u  •  7j,  i  s  Xr  II  7u„  II  , ( w 1 7o ,  I  Z  *  -  I  3  ,  I  Z  )  and  choose  u  such  that 
do  I  Z  0  »  ,.r  1  -j  c 


1  ‘T  e  ^  t|7uo!l»,qJ  *  v  Moreover  H a2 II ‘2 ,3  s  C  J  l7a2lZdxdy  holds 

_U1 

and  we  also  require  1  -Alcs4e  ^  l|7u  II  „/(min  1 0 1  ) t  The  smallness  of  r 

Z  0  <•  _  » 

(see  (4.16))  and  the  smoothness  assumption  on  Tuq  imply  that  there  is  a  y  >  0  such 
that  both  inequalities  hold.  By  proceeding  analogously  with  (in  the  last 

integral  of  l  (3p32>  3j,j2))  the  mixed  terms  are  estimated  below  s  :cn  that  L  is 
coercive  wi th a coerc i v i ty  constant  independent  of  4  and  iU-/Ujl.  Therefore 


(4-22)  11(0,  c  F„(  0  ,0 )  )  '  lI| .  s  const 

Lu ’  v  1,5 

(°c  r  denotes  the  Freche t -der i va t ive  with  respect  to  (E  ,E  )).  Since  Dr  ,  F 

Li*V  U  V  t  »  C  0 

is  uniformly  L i p s ch i tz co n t 1 nu ou s  the  implicit  function  theorem  u  7 

implies  the  existence  of  a  locally  unique  solution  E  u  ,  E  v  of  F5(EU,EV)  *  0  for  s 
sufficiently  small  and 

«(EU’V»1,5  * 

holds.  Then 

(4.23)  'I Eu ::i , 2 ,r.  *  11  £v  "l , 2,r.  3  0{<) 

is  easily  established  Dy  considering  the  probl  ems  ( 4  .  1 9  )  ,  (4.20)  in  qQ  and 
separately. 

To  obtain  an  l<”-es  t  i ma  t e  we  write  (4.19)  as 

div  (aCE^)  »  divf 


Oe  *  in 


4  VJ1  l  Z  ■ 

'  exD(— rrrr  *  ‘  ,n  -;o 
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0 


in  !3 


I  1 

1  .  U  -U,  ,  H 

54exp{-^-r-1 )  5  (•  M)7ufl  in  &c 


(note  that  we  used  ( 4 . 1 2 ) ( a )  to  obtain  f  in  this  form).  The  maximum  principle  in  the 
form  given  by  Gilbarg  and  Trudinger  (1977,  Chapter  8,  paragraph  5)  yields  (since 
f  £L-(n)  and  E,  I  ^  -  0) 

"Eu"-.fl  «  C(a.a)ii7u0i.iao  •'•E'-«'q(ao 


for  any  q>2.  The  embedding  theorem  gives 


II  £ 


n<  const -iwa0H,.a 


"W'l.z,** 


and  (4.23)  implies  ( 4 . 1 7 ) ( a ) .  ( 4 . 1 7 ) ( b )  follows  analogously  and  ( 4 . 1 7 ) ( c )  is  obtained 
from  (4.18).° 

The  implicit  function  theorem  could  not  be  applied  if  there  is  an  n  or  p  region 
without  a  contact  since  uQ  or  vfl  resp.  is  not  uniquely  defined  in  this  region.  The 
linearized  problem  has  a  zero  eigenvalue  then. 


(3.31)  and  (4.17)  imply 

3n  -  0(7) ,  *  °(3> 

The  reduced  current  densities  are  small  close  to  thermal  equilibrium. 

Theorem  4.2  can  easily  be  extended  to  the  case  that  two  n(or  p)  regions  have  a  joint 
boundary  .  The  limiting  function  for  *  (see  (4.17)(c))  remains  unchanged. 


8 )  The  Internal  layer  Problem 

We  now  investigate  the  layer  problem  (3.5)  at  a  fixed  curve  r  of  discontinuity  of  0 
and  prove: 

Theorem  4.3  .  Let  the  assumptions  on  D  on  il  given  in  Theorem  4.1  hold.  Then  (for  given 
i.u.v)  the  internal  layer  problem  (3.5)  has  a  unique  piecewise  monotone 
solution  jiforevery  sEr  which  fulfills  for  every  0  <  *>  <  1 

(  4  .  26 )  (  a )  li(T,s)l  <  C^expf  (  l-«)  •  (0- ^)tor(0-5)t  ♦  0,/  I  ( 0+  ,  s  )  -  »  \  0- .  s  )  I  ) 

for  t  <  -E  *'lj‘  (0*,s)  -  t>‘  (0-,s)l  and 
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(4.26)  (6)  iJ(t  ,s)l  *  C^exp^  (  -1+ui )  /(ir  (  0+  ,"s)  ♦pr  ( 0+  .  s  )  r  +  Q./\  *r(  0+  ,s  ) -f  r  ( 0- ,  s  ;l] 

/Tr 


for  t<E  *  U  (0+,s)-u>  ( 0  -  ,  s )  I  where 

u> 

Cu,0u,Eu>0  depend  on  w  Out  not  on  [  *]  r .  Also  #6C((-*,0)xr)nc((0,«)xr)holds. 


Piecewise  monotone  means  monotone  on  (-*,0)  and  on  (Q,«). 


We  do  not  give  a  proof  for  Theorem  4.3  since  It  follows  completely  the  lines  of  the 
proof  Of  Theorem  3.3  in  Markowich,  Ringhofer,  Sel berherr and  Langer  (1982  b).  We  only 
remark  that 


(4.27)  J(0*.s)  -  O'  (0-.s)(»r(Qr,s)-»~‘  (0-,s))»  n‘  (0  -,  s  )-n  ( 0+  ,  s )  ♦  p1  (0-,s)-pl  (Qf.,sJ 

0r (0+,s)  -  0r(0-,s) 


-r , 


=  r( 


(4.28)  J(o-, S)  .  0r(0>.t)('vr(0*.s)-Jr(0-,s))^  nf  0- ,  s  ) -n1;  0+  ,  s  )  »  $r  (  0- ,  s  ) -pr  ( 0*  ,  s  ) 

0‘ (0*,s)  -  O' (O-.s) 


holds. 


The  estimates  (4.26)  can  be  used  to  compute  the  width  d‘(s)  *  d  +  (  s)+d‘_(s)  of  the 
internal  layer  of  *  at  s  £  r 


(4.29)(a)  d ; ( s )  -  0 


(4.29)(b)  d’.(A)  »  0 


...  * - J 

1  n 

X 

^>n‘  (0+ , S ) rpr(0+  ,  S  )  l 

.  /n  r(0+  ,S  )+p‘  (0+  ,  s  ) . 

[  K  . ( 

1  n 

^r(0-,s)+pr(0-,s)l 

.  -  nr  (0-  ,s)+p'  (.O-.s). 

♦  v7l  i  r  ( 0  +  ,  S  )  - 1 r  ( 0  -,s )  I 


/T7 r 


A  *  ( 0+ ,  s )  -♦  ( 0-  ,  s  )l 


where  d^(s)  (d[(s))  is  the  portion  of  the  layer  width  which  is  on  that  side  of  r  for 
which  t>0  ( t  <  0 )  holds.  (4.29)  is  uniform  in  li).  .  The  width  of  the  depletion  layer 
depends  linearly  on  the  square  root  of  the  potential  drop  accross  the  junction. 


d.,d  can  be  estimated  by  using  the  a-priori  estimates  on  n,  p,  u  and  v.  We  obtain 

IVmaxI 


(4.30) 


101  «  n  *  p  s  101  ♦  0(J*e  . 


If  sgnO‘(0+,S)  «  -sgnOXO-.s)  (that  means  r  is  a  pn-junction)  then  (4.6)  yields 

r , 

1  I  U 

4.31) 


IQ1 (0*,s) -0‘ IQ-.s)  I 

(  0  ♦  ,  S )  -*(0-,S)  «  In  | 

4  IV  1J 

max 

-  In  v'(  0  ,  S  )  +  0  (  5 d  e  ^  ). 


-In  u  ( 0 , s  ) 


If  sgn  0'(0*,s)  *  sgn  O'  ( 0 - , s )  ( r  is  an  nn  or  pp  junction)  then  (u]„  is  uniformly 


bounded  as  5  ■*  0. 

(4.32)  i[(t)  • 


We  get  from  (4.29) 


(  ■  x  f 

,  A 

1  n  - 

\/0r(0+,s)  V 

/Dr(0+,s) 

■W 

4 

▼  o  6 


(if  r  is  a  pn-junction) .  An  analogous  formula  holds  for  d[(s),  The  square  root 
drops  out  if  r  is  an  nn  (or  pp)  junction. 


I  If 


maxi  \ 

*r 


C )  The  Semi conductor-Ox i de  Interface  Problem 


Theorem  4.4  The  problem  (3.20),  (3.25)  has  a  unique  monotone  solution  *  (for  every 

q  €  3a..)  which  fulfills  for  every  0  <  u  <  1 
os 

(4.33)  l*(o,q)l  «  Cuexp((-l  +  w)v,n3-(0.q)+p5(0.q)o+  0u/)?(0,q)l 

for  o  >  E  /l  *(  0  ,q )  I  where  C.,,D  ,£  >0 
are  independent  of  7(0, q).  Also 

(4.34)  U(0,q)l  <  n(  U3(0,q)  l2,q) ,  (b)  sgn*(0,q)  *  -  sgn7p(0,q) 

holds,  where  n(-,q):[0,«)  -  [0,«)  is  a  monotonely  increasing  function  which 
fulfills 

n(t  ,q )  ~  -  y* - - —  JS  t -0  and  n(».q)  *  «*,  nt(<-,q)  *  0. 

/i3{0.q)+p3(0.qT 


Proof :  From  Fife  (1974)  we  derive  that  there  is  a  unique  monotone  solution  of  (3.20) 
with  given  boundary  value  7(0, q).  At  first  we  assume  that  <i(0,q)  >  0. 

Then  *( a ,q)  fulfills 


( 4 . 35) ( a )  a 


*{0,q) 

J 

*( o ,q) 


_ d  a 

✓2S(“  Jq) 


where 

(4.35)  (b)  S(  i  ,q )  •  /(S^O.qle1  -  p^O.qJe*1  -  D3  (  0  ,q  ;  )dt 
0 

.  n  5 ( 0  ,  q ) (ea*l )  ♦  p3(0.q)(e'a-l)  -  03(0,q ) a 


holds  (see  Fife  (1974)). 


Differentiating  (4.35)(a)  with  respect  to  a  and  evaluating  at  a* 0  gives  the 
equat i on 

(4.36)  -*  fO.q)  *  v'2G(*(0.q),q). 

P 

Therefore  the  interface  problem  (3.20),  (3.25)  has  a  unique  monotone  solution 
for  e  (0,q)  <  0  if  (4.36)  has  a  unique  positive  solution  *(0,q). 

For  0,q)  <  0  we  apply  the  same  argument  to  -1  and  find  that  (3.20),  (3.25) 
has  a  unique  monotone  solution  for  “JT  (0,q)  >0  if  the  equation 

(4.37)  v  (0,q)  »  /2G(«(0.q) ,q) 

J 

has  a  unique  negative  solution  uT ( 0 , q )  .  From  *  (0,q)  *  0  we  immediately  get 
»(u.q)  *0.  A  simple  calculation  shows  that  the  function 

H ( a ,q )  =  i 23( a  ,q  ) 

is  decreasing  for  j£(-»,01,  increasing  for  a  E  [0,»)  and  H(-«,q)  «  H(«,q)  *  =» 
for  every  q  6  3 P2Q s  . 

This  settles  the  existence  statement  of  the  theorem  and  (4.34)(b). 

The  estimate  (4.33)  follows  by  proceeding  as  in  Markowich,  Ringhofer,  Selberherr 

and  Larger  (  1982  b.  Theorem  3.3).  To  prove  (4.34)  ( a )  we  set 
2 

l  *  iLi — { and  rewrite  (  4.36),  (  4.37  )  as 
2 

G( Z  ,q)  *  3  ,  Z : «  *(0,q) 

for  fixed  q.  At  first  we  investigate  the  root  Z^  >0.  D i f f erent i a t i on  with 
respect  to  a  yields 

jf  (Z1(3 )  ,0)Z{(S  )  *  1  . 

Zj(®)  6  ( 0 ,«)  implies  Zj(°°)  €  (0,a>),  which  is  a  contradiction.  Therefore  Z^(®) 
and  -f^(<=,q)  =  ^implies  Zj>l  =  0.  The  root  Z,  <0  is  investigated  analogously 
(Zjf'”)5-®,  Zj(«)»0  hold). 

Taylor  expansion  of  5  for  2  -0  gives  the  behaviour  of  Z  close  to  zero.o 

The  width  d^Jq)  of  the  interface  layer  is  . . .  •  ^ 

'max 

_ __  T  i 

(4.33)  dg$(q)  «  0.  -  ■  ■  —  !  I  n  —  ‘ -  *  v'l7(0,q/l  *  S^e  T 

<  v7!  0  5  (  o  t  a )  ■  i i  |  /|  o  5  ( 0  ,  q )  i 
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The  first  order  term  of  the  expansion  «  +  $  +  1  cannot  be  expected  to  approximate 
the  solution  of  the  singularly  perturbed  problem  close  to  the  point  S  with 
{ S  >  *  F  n?nQs  since  1(0, •)  is  discontinuous  at  q  »  S  unless  7(0,-)  •  0. 

Generally  the  first  order  term  can  only  approximate  the  solution  outside  a  small 
sphere  (with  radius  at  least  0(x(1nx())  around  S. 

Theorem  4.4  can  be  used  to  compute  an  estimate  for  the  threshold-gate-voltage,  thatisthe 
gate-voltage  at  which  strong  inversion  (the  minority  carrier  density  at  the  inter¬ 
face  is  at  least  as  large  as  the  absolute  value  of  the  dooing  profile)  starts. 

To  demonstrate  this  we  take  the  n-cnannel  MOS  transistor  as  shown  in  Figure  1.  n  is 
the  minority  carrier  density  in  the  p-region  therefore  strong  inversion  occurs  if 

n3(0,q)  »  1 0 3 ( 0  ,  q ) I 

holds  for  some  q  £  JG  .  from  ri>  '  j  we  (up  t0  O(X)-terms)  n3(0,q)  « 

««  n3  (O.q)e11^  such  that  the  strong  inversion  condition  is 


1(0, q)  >  In 


i  p  ■;  o ,  g ;  i  ^ 
n  3  (  0  ,  q  )  1 


~tnv/<°  •<!) 


(4.36),  (4.37)  yields 


,Vmax’ 


4  — rjj 

((1,(0, q)  !2  »  2G(1,NV(0,q),q)  s  2  - - 

3  I D3 ( 0 , q )  t 


-2l03(0,q}lln  +  2  I  D 3  (  0  ,  q )  I- — 

6  T 

and  from  (3.26),  f  2 . 8 )  ( b )  we  derive  that  strong  inversion  occurs  if 


UG  >  UF  +  des 

*TT0 


« -  v-r 

5  e  T 


max  I  0  I 

a, 


♦  2max 1 0  l  1 n 


+  2max I D I 


"'max' 


holds. 


( 


5.  ASYMPTOTIC  REPRESENTATION  OF  SOLUTIONS 

To  demonstrate  the  validity  of  the  expansions  we  show  that  if  the  device  is  In 
thermal  equilibrium,  then  the  solution  of  the  singularly  perturbed  problem  is  close 
(in  l"(n))  to  the  zeroth  order  term  in  the  asymptotic  expansion  (3.1)  assuming  that 
x  is  small.  For  simplicity  we  assume  that  there  is  only  one  exjunction  r 
( a  •  n0  u  r  U  cjj) .?  intersects  with  an  in  two  points  S^.Sj  which  have  positive  distance 
from  Also  there  aresphere.s  with  radius  a  centered  at  and  S2  within  which  r 

and  3  a  are  perpendicular  lines. an  is  assumed  to  be  as  in  Theorem  2.2.  Moreover  no 
semiconductor-oxide  interface  occurs,  i.e.  «  *(  )  . 

We  have  to  overcome  one  technical  difficulty,  namely  that  the  layer  solution  #(•£, s ) 
is  only  defined  in  an  open  u-strip  S  (r)  about  r.  In  order  to  extend  this  function  to 

U 

a  we  choose  $eC“(R)  with 

8  ( t )  *  1  for  t  €  ( -  2  ■  2  3  *  for  ltl*u 

and  set 

S(!L*zll,  S(x,yl)  •  3(t(x,y))  ,  (x,y)  £  SJr.) 

(5.1)  S(x,y))  J 

1  0  ,  (x.y)  £ S  (:) 

k  a 

(see  Fife  (1974). 

We  prove: 


Theorem  5. 1  Let  r  and  a  fulfill  the  aoove  stated  assumptions.  Also  D  is  piecewise 
constant ,  that  means 

r  0  in  a 

D  *  i  „  where  D„,D,  are  constants  with  D„  *0,.  Then 

i  u  1  1  n  u*  01  0  1 

1  1  3 

(5.2)  *(x,y,A)  «  u(x,y)  +  »  (  ,  s(x,y))+  0  ( Alina  1^). 

e  * 

5 

7 

holds  if  x,  xilnfi  are  sufficiently  small  and  if  the  device  is  in 
thermal  equilibrium  (i.e.  U^«0) 

Proof :  In  thermal  equilibrium  u  •  v  ■  1  holds  and  (2.4)  simplifies  to 
2  2 

(5. 3)  (a)  k  a*  ■  Zs  -sinh*  -0  in  .0 


subject 

I5.3)(b) 

We  set 


to  the  boundary  conditions 


3r  "  3. 

"is  "C 


,.!  O  +  v'd2*!,4  1 

1  r, - , - 

v  1 
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(5.4) 


<i»  *  Y  ♦  <0 


where  the  reduced  solution  I  is  given  by  (4.6) 


(5.5)  J  •  1 n 


0+v/02+4«4 


since  u«v«l  holds.  *e  is  given  by  (5.1)  (*(t,s)  solves  (3.5)). 

Assuming  that  u  is  less  than  the  distance  of  Sj.Sj  from  an^.  we  get 
(  5 . 6 )  (a )  «|Jfl  •  0 
and  because  of  (5.1) 

( 5 . 6)  ( b)  ?<p-«l3n  •  0 

(note  that  7*-n!  *  0  because  0  is  piecewise  constant  and  r  is  perpend i cu 1 ar  to 

an).  s 

Inserting  (5.4)  into  ( 5 . 3 ) ( a )  gives 


I  T  ♦  a  J  it  +  V^dio  *  26“  -  si  nh  (i+i+u>)-0  in  S  (')  ■ 

2 


(note  that  i(J*J  )  (  i  n  the  weak  sense)  is  in  L*('.)  and  equals  -2»eTT+  s’nce  *•» 

do  not  depend  on  Sand  since  the  interface  conditions  (3.5)(cj,(d)  hold).  Also 


holds. 


1  A  n  n  •  /V 

\  i*„  ♦  \  iip  »  2S.sinh(**w*<p)-0  in  S  (r)  -  S  (r) 

e  e  a) 


2  2 

x  dip  *  2a  "s i nh( *  *<o)-D  m  n-s  (r) 


By  using  the  internal  layer  equations  (3.5)(a),(b)  we  rewrite  (5.7)  as 

-  -  A 

(5.10)  \  iip  J  25  cosh(ji+*)qi  +  J  e  *+  **  (  e^-<p- 1 ) 

.  42g-I-*(  )  -  x*  it  in  S  (r) 


and  by  using  (3.2)(a)  (with  u«v»l)  we  obtain  f-om  (5.8),  (5.9) 


I 


,5. Ill 


.‘"ios  *  252  cosh(**®  )<p  +  iZe“  i'e  ( e^-g- 1 ) 


■>  '  A 

-  5£e'*"'e (  e"'®+<p- 1 ) 


-  iV^e  *e-l)  »  62e*(e’’e-l)  -  v2dj.  in  S,(r)  -S 


(5.12)  \Zm  *  252cosh(  J  )g  ♦  42e*(e'*-g-l)  -  4  2e‘*(  e'^+g- 1 )  in 

Therefore  we  have  a  oroblem  of  the  form 

(5.r1)(a)  k2.ag  »  2  5  2cosh(i+C  )g  ♦  F ( g,x ,y , x , a )  in  a 


(5.13)(b)  £  I  *  7g- *  0 

3“C  ’-is 


wnere  F  fulfills  the  estimate 


,  2 


(5.14)  !  F  ( 0,  x  ,  y  ,  l ,  ? )  i  s  C(k)(g  *llt'e!laaS  (r)-S  s  , 

‘  Ji  /  -  v  k  *  vJ 

jf.  )  2 


''^V1  »,s  s  ), 


(*,y)  e 


*"or  igisk  with  C .  < )  independent  of  5  anO  v. 

F'om  Markowich,  Ringhofer,  Selbernerr  and  langer  (1982b)  we  get 

I  a.  5  s  maxu  <•  .(0*n  .  U  T(0  -  )  I )  s  D,  »  1 1  ns  I 

2»' :) 

and  ( 4 . 26)  ,  : 5.  1 )  imply 


I  «jl 


r  s  Oy— T  ®  x  P ( ~  T ) 


e"-,S  !D-S  ID  s  U27T 
"  j  5 

(D-S  ID  *°3~r7  exp('"J 

J  *  2 

2 

wnere  3,  ,02.03,2  >  0  are  independent  of  i  and  5  as  \-Q+,  i-  0*. 
This  yields 


0  .  :  4  )  1  F  (  g ,  x  ,  y  ,  V  ,  J  )  !  s  C ,  .  k  )  !  g2  *•  1  .  "iT'  4  I ) 

for  igisk,  i,5and  •  1 1 n  5  1  sufficiently  small. 

We  now  drove  in  (."-estimate  for  the  solution  of  the  oroblem 


(r) 

n  -  S 
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( 5 . 15) (a)  Lx>Jw  « 

(  5  - 1 5  )  (  b )  w  ■  *  Vwn  I  ,  *0 

5"C  :,!is 


2  7  -m  A 

rj«  -  2s  ‘cosn{  it»*  )w  ■»  f(x,y)  in  .i 


with  f  e  i*(a). 


Therefore  we  define  the  function 

(5.16) 

with  L  >  0 ,  3  >  0  • 


Vb(x.y) 


+  exp  — p — -  t2(x,y)|  9(t(x,y)) 

\  x  c  I  1  n  6  : 


A  straight-forward  calculation  (analogous  to  the  proc-f  of  Theorem  4.3  in  Markowich, 
Ringhofer,  Se 1 berherr a  no  Linger  (  1982b) )  which  heavily  relies  on  the  estimate  (4.26)  show:, 
that  the  constants  1,3  can  be  chosen  i ndepenoen t 1 y  of  such  that 

s  -  TtSIt  •  !*•*>“ 


»  L  , 


'V<n. 

is 


holds.  This  implies  that  V &  can  be  used  as  compansm  ‘unction  and  the  maximum 
principle  yields 


(5.17)  I  w||  «  F  1 1  nS  |  !|  f  II  , 

oot tip,  .4 

for  the  solution  w  of  (5.151  (  F  >  0  is  independent  of  1,5). 

3 

7 

7  2  >  ^ 

Now  we  define  tne  operator  ?* :  A  _  .-L'(r.)  where  ^  *  i  so  €  U  ( it )  I  1 1  ip  II  ^  s  5xi'.n5l 

with  G  >0  independent  of  i  and  5  such  that  so= M  (  3 )  is  the  solution  of 

0  -X 

t'dip  3  2  Cosh(  ***  )<®  ♦  F  (  0  ,  s  ,y  ,  v ,  5)  in  n 

U)  *  dip  ’  rt  *  0  . 

1  -‘o  H  s 

A  ‘ixeo-oointio*  of  M  is  a  weak  solution  of  (5.13),  since  Range  (H)  crl  ‘(.5).  M  is 
continuous  and  therefore  Ascoli's  theorem  implies  that  it  is  compact.  Fnom  the 
estimates  (5.14),  (5-17)  we  conclude  that  M  map  A x  ,  into  itself  if  G>C|{ k)F  and 

,  •  1  n  5  1  ^  i  s  suf  f  i  C 1  ent  1  y  small. Since  A  is  closed  3nd  convex  in  L  ‘( )  Schauder's  fixed 
point  theorem  yields  the  existence  o'  a  solution  of  (5.13)  and  (5.2)  follows.  As 
mentioned  before  this  solution  ;s  unique  (in  h  ^ n )  )  .  o 


~  neorem  5.1  can  easily  be  extended  to  doo • ng  profiles  wnich  are  not  piecewise  constant, 
'he  proof  proceeds  aiong  tne  same  lines  eicept  t.nat  A( **»e )  -  1  °°( .1 )  if  ■>  •  ve  depend 
on  s.  'his  can  be  overcome  by  modifying  tne  interface  condition  (3.5)(d) 
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I 


appropri ately ( see  Markowich,  Ringhofer,  Selberherr  and  banger  (1982b)  for  the 
one-dimensional  problem).  The  extension  to  more  than  one  junction  is  also  straight 
forwa rd . 


A  main  assumption  of  the  Theorem  is  that  xilnl 


sma  1 1  . 


(2.6) 


implies  that 


(5.18)  s2  ■  y 2 x 2 ,  y2  • 

holds.  For  realistic  devices 
proportional  to  x  and  xilnil 


5  y  *0.1.  therefore  6  can  be  regarded  to  be  at  least 

7  ^ 

~  xllnx|2  as  x  - 


0+. 


6.  NUMERICAL  RESULTS  AND  EXTENSIONS 


We  present  numerical  results  for  a  two-dimensional  diode.  The  geometry  of  the  device 
is  depicted  in  Figure  2. 


Figure  2.  Oiode 


(2.6)  gives  ^  *  4^  *  10’7 . 

The  computations  described  in  the  sequel  were  performed  using  a  specialised  finite 
difference  d i scret i sa t ion  method  completely  based  on  the  singular  perturbation 
approach  (see  A.Frant  et  al  (1982)). 

Figure  3  shows  the  potential  in  thermal  equilibrium  (UQ*Uj*0V),  Figure  4  shows  the 
electron  density  n  and  Figure  5  the  hole  density  p.  As  proven,  in  Section  4  the 
n-region  is  (almost)  depleted  of  holes  and  the  p-region  ,Tg  is  (almost)  depleted 
of  electrons.  The  internal  layer  at  the  pn-junction  fj  is  clearly  visible.  Figure  6 
shows  the  grid  used  for  the  d 1 scret i sa t ion  (which  was  generated  automatically  by 
the  code).  An  accumulation  of  grid  points  around  the  pn-junction  (in  order  to 
'balance'  the  large  derivatives  of  the  solutions  in  perpend icul ar  direction  to  Tj) 
is  evident. 

Figure  7  shows  the  potential  distribution  for  -10 V  applied  bias  (i.e.  U*UQ -U ^ *  -  10V ) . 
The  width  of  the  depletion  layer  increased  (compared  to  the  equilibrium  problem) 

(see  (4.32)). 

The  electron  density  n  for  U«0.75V  applied  bias  is  depicted  in  Figure  8. 

4  - 1 

5  exp( I Ul / UT )  « 10  (we  are  dealing  with  a  low  injection  case)  and  n  visibly  'lifts 

-  4 

off’  from  zero  in  the  p  region  (the  theory  predicts  n » 0  ( 5  exp(  I U I /Uy.) )  inro). 

4  3 

Figure  9  shows  the  hole  density  p  for  U*1V  applied  bias.  Now  6  exp( I U 1 /U y)  «« 2 .  3 5  x  10 
holds  ( h igh- i n jec t i on ) .  p>>IDIaway  from  the  contact  in  the  p-region  <ig  and  the 
n- region  is  also  flooded  with  holes.  The  absolute  value  of  J?,  that  is 
I J p I  »  i/( J* ) ^  +  ( J y ) ^  ,  is  depicted  in  Figure  10.  The  'peak'  at  the  edge  of  the  contact 
C0  represents  the  singularity  at  intersections  of  J.Tj.  and  ”•  1S  as  discussed  in 
Theorems  2.1,  4.1.  This  phenomenon  is  physically  interpreted  as  current  crowding  at 

contact  edges. 

» 
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Figure  1C 
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In  this  paper  we  only  dealt  with  abrupt  doping  profiles  {that  means  D  has  discontinuities 
accross  junctions).  This  assumption  is  actually  a  simplification  since  'physical' 
doping  profiles  are  exponentially  graded  near  junctions.  The  singular  perturbation 
approach  can  be  modified  to  accomodate  this  more  realistic  situation  by  setting 

(6.1)  0(x,y,x)  •  B(x.y)  ♦  S( » .  s(x,y)) 

A 

{close  to  a  particular  junction  r)  where  5  is  abrupt  and  ID(r  .  s ) 1  <  £ ( s ) exp( -F ( s ) t ) 
holds  with  £,F>0.  0(.,.,\)  is  continuous  along  r  for  \  >  0. 

A 

Then  0  appears  as  'inhomogeneity'  to  the  internal  layer  equation  ( 3 . S) ( a )  ,  ( b) .  Other¬ 
wise  the  preceding  theory  remains  unchanged  when  0  is  substituted  for  D. 

As  mentioned  in  the  introduction  we  neglected  recombination-generation  of  carriers  in 
our  model -equations  (1.1).  The  recomb i na t i on -genera t i on  rate  R  in  (l.l)(b),(c)  is 
(in  the  most  general  setting)  a  nonlinear  function  of  n,p,Jn,Jp  and  do  (see  Sze  (1981) 
for  details). 

The  existence  proofs  (Theorems  2.1  and  4.1)  are  affected  by  the  introduction  of 
recombination  (since  they  use  the  linearity  of  the  continuity  equations  in  u  and  v). 

The  qualitative  behaviour  of  solutions  is  pretty  much  unaffected  in  1 ow- i n j ec t i on  con¬ 
ditions  since  R«0  in  thermal  equilibrium.  In  hi gh- i nj ecti bn  conditions  however 
recombination  generation  has  a  decisive  impact  on  the  solutions  and  cannot  be  neglected 
(see  Schiitz,  Selberherr  and  Potzl  (1982)  for  a  numerical  study  of  'avalanche'  effects 
i n  MOS- t  rans i S  tors ) . 

Tne  validity  proof  for  the  asymptotic  expansions  in  non-equilibrium  is  an  unsettled 
issue.  For  the  one-dimensional  semiconductor  problem  Markowich,  Ringhofer,  Selberherr 
and  Langer  (1982  b)  gave  a  proof  (with  estimate  of  the  remainder  term)  for  a  diode 
close  to  thermal  equilibrium.  There  is  numerical  evidence  that  the  asymptotic 
expansions  'represent'  a  solution  even  for  large  applied  voltages  (see  Markowich, 
Ringnofer,  Selberherr  and  Langer  (1982a, b)),  but  no  estimate  of  the  remainder  term 
for  arbitrary  bias  is  known  so  far. 
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